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Large scale structure introduces two different kinds of errors in the luminosity distance estimates 
from standardizable candles such as supernovae la (SNe) - a Poissonian scatter for each SN and a 
coherent component due to correlated fluctuations between different SNe. Increasing the number of 
SNe helps reduce the first type of error but not the second. The coherent component has been largely 
ignored in forecasts of dark energy parameter estimation from upcoming SN surveys. For instance 
it is commonly thought, based on Poissonian considerations, that peculiar motion is unimportant, 
even for a low redshift SN survey such as the Nearby Supernova Factory (SNfactory; z = 0.03 — 0.08), 
which provides a useful anchor for future high redshift surveys by determining the SN zero-point. 
We show that ignoring coherent peculiar motion leads to an underestimate of the zero-point error 
by about a factor of 2, despite the fact that SNfactory covers almost half of the sky. More generally, 
there are four types of fluctuations: peculiar motion, gravitational leasing, gravitational redshift 
and what is akin to the integrated Sachs-Wolfe effect. Peculiar motion and leasing dominates at low 
and high redshifts respectively. Taking into account all significant luminosity distance fluctuations 
due to large scale structure leads to a degradation of up to 60% in the determination of the dark 
energy equation of state from upcoming high redshift SN surveys, when used in conjunction with 
a low redshift anchor such as the SNfactory. The most relevant fluctuations are the coherent ones 
due to peculiar motion and the Poissonian ones due to lensing, with peculiar motion playing the 
dominant role. We also discuss to what extent the noise here can be viewed as a useful signal, and 
whether corrections can be made to reduce the degradation. 

PACS numbers: 98.80.-k; 98.80.Es; 98.80.Jk; 95.30.Sf 


I. INTRODUCTION 

The problem of the cosmological constant, or more gen¬ 
erally dark energy, is one of the deepest problems in cos¬ 
mology today. While there are by now multiple lines of 
evidence for the existence of dark energy Q , the evidence 
from type la supernovae (SNe) was historically what con¬ 
vinced a large fraction of the cosmology community that 
this enigmatic form of energy should be taken seriously 
iQ Upcoming and ongoing SN surveys 0, with vastly 
improved statistics, promise to constrain the equation 
of state of dark energy to unprecedented precision, thus 
shedding light on the issue of whether the apparent ac¬ 
celeration of the universe is caused by the cosmological 
constant, a dynamical scalar field or departure from Ein¬ 
stein gravity Q- 

There has been much recent work on projections for 
the determination of dark energy properties from these 
SN surveys. By and large, they focus on the following as¬ 
pects of the error budget: intrinsic statistical error, sys¬ 
tematic error and grav itational lensing induced scatter 
(e.g. [^I^I^I^ ITfl illL ll ^ and references therein). The 
intrinsic statistical error refers to the intrinsic spread in 
SN luminosity even after suitable standardizing correc¬ 
tions have been applied. It is typical to assume that the 
intrinsic spread in magnitude has a (root-mean-squared; 
rms) size of cr™*''' = 0.1 — 0.15 for each SN ^3- This 
kind of intrinsic statistical error can be beaten down by 


having a large number of SNe. There are several sources 
of systematic error, such as Malmquist bias, luminosity 
evolution, imperfect corrections for dust extinction, and 
so on. They are not necessarily diminished by having a 
large number of SNe, although a large sample often helps 
in identifying and characterizing them. Lastly, gravita¬ 
tional lensing by intervening structures introduces fluc¬ 
tuations in the observed flux of SNe. So far, the focus 
has been on how gravitational lensing introduces a Pois¬ 
sonian scatter rather analogous to the intrinsic spread. 
This kind of error can likewise be reduced by having a 
large sample of SNe nm. 

The existing discussion can be improved in two ways. 
First of all, gravitational lensing by large scale structure 
introduces not only a Poissonian scatter to the individual 
SN flux, but also correlated flux fluctuations between dif¬ 
ferent SNe. One can view the correlated fluctuations as a 
consequence of the large scale coherence of the interven¬ 
ing structures. Second, large scale structure introduces 
fluctuations beyond that captured by gravitational lens¬ 
ing, and like lensing, these fluctuations have a Poissonian 
component as well as a correlated or coherent component. 
It is worth noting that an expression for all the hrst or¬ 
der fluctuations in the luminosity distance - hrst order in 
metric and energy-momentum perturbations - has been 
worked out for quite some time e.g. miniii (with 
minor corrections; see below). The full implications for 
current and future SN surveys, however, have not been 
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explored, with an important exception (pointed out to 
us by Dragan Huterer) - Sugiura, Sugiyama & Sasaki 
M computed the anisotropies (the dipole and beyond) 
in luminosity distance and investigated the implications 
for measurements of the decceleration parameter go- 

As we will see, to first order, there are four sources of 
luminosity distance (or magnitude) fluctuations: gravi¬ 
tational lensing, gravitational redshift, peculiar motion 
and an effect akin to the integrated Sachs-Wolfe (ISW) 
effect. We will see that for most practical purposes, it 
is sufficient to consider gravitational lensing and peculiar 
motion. They become important at high (z > 1) and low 
{z < 0.1) redshifts respectively. 

What is particularly interesting, and perhaps surpris¬ 
ing, is that peculiar motion plays a significant role in the 
degradation of dark energy errors. There is a widespread 
perception that the effects of peculiar motion are neg¬ 
ligible as long as the median redshift is greater than 
0.05 or so. Let us take as an example the Nearby Su¬ 
pernova Factory (SNfactory; other low redshift surveys 
include the CfA Supernova Pr^ram, Carnegie Super¬ 
nova Project and LOTOSS, see 0), whose redshift range 
z = 0.03 — 0.08 was chosen in the hope of making the ef¬ 
fects of peculiar motion negligible. Such a perception 
seems at first sight quite reasonable: typical peculiar ve¬ 
locities are of the order of 300 km/s, and so the ratio 
of peculiar flow to Hubble flow at z = 0.055 is about 
300/(3 X 10® X 0.055) ~ 0.02. Translating this into fluctu¬ 
ations in magnitude (details are given in SZZI), we have 
6m ~ 2.17 X 0.02 ~ 0.04, which is quite a bit smaller 
than the intrinsic spread in SN magnitude (0.1 — 0.15), 
apparently suggesting we can ignore peculiar motion (re¬ 
call that different sources of errors are to be added in 
quadrature). 

What such an argument misses is that coherent pe¬ 
culiar flows introduce correlations in magnitude fluctua¬ 
tions between different SNe. While it is true that peculiar 
motion introduces a negligible Poissonian scatter com¬ 
pared to the intrinsic scatter, the correlated component 
cannot be ignored as it turns out. One can intuitively 
understand it as follows. As the number of SNe {N) be¬ 
comes large, the intrinsic statistical error is beaten down 
to be quite small in the usual root-N fashion. Correlated 
errors, such as that due to correlated/coherent peculiar 
flows, are not reduced by N at all, and so there must 
be some N beyond which the correlated errors become 
dominant. We will see that this is indeed the case for 
the SNfactory. Since a low redshift survey such as the 
SNfactory plays an important role in constraining the SN 
zero-point, dark energy determination from higher red¬ 
shift surveys (where peculiar motion is less of an issue) 
is affected indirectly by these considerations as well. 

Coherent large scale flows (i.e. bulk flows) have of 
course been the subject of research for a long time (see 
M for a review). Particularly relevant to our investiga¬ 
tion are discussions of the peculiar velocity monopole, or 
what is sometimes referred to as the local Hubble bubble, 
which incidentally made use of SNe la ESmiH. We 


will see later that for a survey like the SNfactory, which 
covers roughly half of the sky, fluctuations in the lower 
velocity multipoles (the monopole, dipole, etc) contribute 
significantly to the dark energy error budget. 

The rest of the paper is organized as follows. In M we 
set the stage by describing the parameters of interest and 
how the fluctuations in magnitude are related to the pa¬ 
rameter errors. Some details on how to go from the Fisher 
matrix to actual errorbars of various types are given in 
Appendix^ We describe in m how the average mag¬ 
nitude in a given redshift bin fluctuates, and how these 
fluctuations can be divided into a Poissonian component 
and a correlated/coherent component, in effect defining 
the magnitude covariance matrix. To keep the discussion 
simple, the derivation is relegated to Appendix m In 
we derive an expression for the first order luminos¬ 
ity distance fluctuations (SZH), and work out explicitly 
their implications for the magnitude covariance matrix in 
terms of the mass power spectrum To keep the 

discussion simple in the main body, details of these two 
steps are relegated to Appendices O and IDI Appendix O 
might be interesting to the more theoretically inclined: 
the explicit expression given here for the luminosity dis¬ 
tance fluctuation corrects a minor error in earlier expres¬ 
sions in the literature. Appendix IHI contains expressions 
for the velocity window function for an arbitrary survey 
geometry that might be of interest to observers who are 
interested in making predictions for their own surveys. 
In gV| we finally put everything together to make error 
forecasts. It proves illuminating to first focus separately 
on the contributions to errors from peculiar motion and 
lensing (issi and EH), and in we make forecasts 
for a number of ongoing/planned/proposed SN surveys. 
The key results are summarized in Fig. [3|H1 Eland Table 
El We conclude in EH with a brief summary of major 
results and a discussion of several issues that naturally 
arise, some of which are worth exploring further: 

• whether peculiar motion degrades the current con¬ 
straints on dark energy (the answer is: not signifi¬ 
cantly) ; 

• how the exact survey geometry impacts dark energy 
errors; 

• whether internal motion that could add to the pe¬ 
culiar velocity is important (the answer is no); 

• the issue of systematic errors, and how they might 
change our conclusions; 

• whether realistic redshift measurements are accu¬ 
rate enough for us to have to worry about peculiar 
velocities (the answer is yes); 

• whether corrections can be made for peculiar mo¬ 
tion and lensing to reduce the dark energy errors 
(the answer is: probably difficult for (the Poisso¬ 
nian part of) lensing, but maybe yes for peculiar 
motion); 
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• whether the noise that we refer to here from pecu¬ 
liar motion and lensing can in fact be turned into 
a useful signal (the answer, for lensing, is that the 
signal is not competitive with the lensing of galax¬ 
ies). 

A comment on our terminology: we refer to the fluc¬ 
tuations of interest in this paper as large scale struc¬ 
ture induced. The term large scale structure should be 
viewed as synonymous with departure from homogene¬ 
ity. Some of the fluctuations discussed here, such as 
the Poissonian lensing scatter, are in fact dominated by 
structures on relatively small scales (galactic scales or 
smaller). Also, even though much of the discussion in 
this paper is phrased in terms of SNe as standard can¬ 
dles, most of our expressions of course apply equally well 
to any other distance indicators. 

While this paper was in preparation, two preprints 
appeared in the electronic archive that partially overlap 
with ours, specifically concerning lensing covariance as 
noise and signal. See also the preprint by on SN 
lensing as a potentially useful signal. There is also a 
preprint by that discusses fluctuations of the lumi¬ 
nosity distance in general as a useful signal (see also [2^1. 


from experiment to experiment, but is generally close to 
but slightly less than unity. Throughout this paper, we 
assume a flat universe, so the matter density rim is not 
an independent parameter. 

To summarize, given a SN experiment yielding m as a 
function of 0 , we can fit for four parameters: iCpivot, Wa, 
ride and M. Let us label them by Pa with a ranging from 
1 to 4. We will mostly marginalize over M since it is not 
of cosmological interests. Sometimes we marginalize over 
ride as well, usually with a prior. Exactly what prior, 
if any, is used will be stated explicited in each worked 
example below. 

The Fisher matrix is defined by 




drrii -i druj 

dpa dpp 


( 3 ) 


where we imagine that the SNe have been binned up in 
redshifts labeled by Latin indices: mi refers to the aver¬ 
aged m for SNe that fall within a redshift bin centered at 
z = Zi- Here a and range from 1 to 4, corresponding 
to the parameters Wpivot, Wa, ^de and M. The (binned) 
magnitude covariance matrix is 


II. PRELIMINARIES 

The relation between luminosity distance (11 and ap¬ 
parent magnitude m is: 

m = 5logj^gdi-l-M = hlndi/lnlO-l-M ~ 2.171n(ii + Af 

( 1 ) 

where M is the magnitude zero-point. Note that it is cus¬ 
tomary to define the absolute magnitude as differing from 
our M by some additive constant, whose precise value de¬ 
pends on the unit used for cIl- If we rescale di by some 
multiplicative factor, such as when we switch units or 
when we alter the Hubble constant today Hq (e.g. if we 
express di in Mpc/h), the change can be absorbed into 
the definition of M. Ultimately, in determining cosmo¬ 
logical parameters from m(z), the apparent magnitude 
as a function of redshift z, we would marginalize over 
the zero-point M, which means marginalizing over the 
absolute magnitude and Hq at the same time. 

What are the cosmological parameters of interest? 
Here, we are interested in fide, WpWot and Wa, which are 
respectively the dark energy density today (normalized 
by critical density), the dark energy equation of state and 
its evolution. We use the parametrization proposed by 

mm: 

w(^Cl) = X/^pivot X/^a(^pivot (2) 

where w{a) is the equation of state of dark energy at 
scale factor a, Wpivot is the equation of state at a = Opivoti 
and Wa is the negative slope —dw/da. The scale factor 
Opivot is chosen such that the errorbars on rCpivot and 
Wa are uncorrelated. The precise Opivot therefore varies 


Cij = (SmiSmj) (4) 

All relevant errorbars related to the four parameters, 
marginalized or otherwise, with or without prior, can 
be deduced from the Fisher matrix Fa/J- For instance, 
the (rms) errorbar on Wpivot marginalized over everything 
else with no prior is given by \/[F“i]ii. Further de¬ 
tails are given in Appendix^ especially on how to choose 
the pivot scale Opivot- One particularly useful fact: as 
long as Opivot is chosen such that the errors on Wpivot and 
Wa are uncorrelated, the error on Wpivot marginalized over 
Wa is exactly the same as the error on w setting dw/da 
to be some fixed value, say zero This is true even if 
there are other parameters present, such as Hde and M. 
To be precise: the error on rcpivot marginalized over Wa 
and other parameters (let us call them pq, pd ...) is ex¬ 
actly the same as the error on a constant w marginalized 
over P 3 , pd ... (see Appendix 0. This fact is useful for 
comparing our results with some of those in the literature 
which often assume a constant w. 

A common alternative parametrization, w = wq + 
w'{z — Zpivot), would result in an errorbar for wq that 
is similar to Wpivot, and an errorbar for w' that is typi¬ 
cally about half of that for Wa- 

Any error projections for surveys necessarily assume 
a fiducial model. We will assume throughout a flat cos¬ 
mological constant dominated model with matter density 
flm = 0.27, dark energy density fide = 0.73, Wpivot = —1 
and Wa = 0 |^ . 

Our next task is to calculate the magnitude covariance 
matrix C^. 
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III. THE MAGNITUDE COVARIANCE 
MATRIX - POISSONIAN AND COHERENT 
COMPONENTS 


Let us recall that we have from observations a vector 
of numbers: the apparent magnitude mi averaged over 
redshift bin with Zi being the average redshift of that 
bin [^. As we show in Appendix m the magnitude 
covariance matrix is 


/ intr.',2_|_C Poiss.')2 

Q, = {6m,Sm,) = ^ + C,, (5) 


where Ni is the number of SNe in the i-th redshift bin, 
^mtr. jg intrinsic dispersion of SN magnitude which 
we will take to be either 0.1 or 0.15 [13, and is the 

Poissonian dispersion induced by large scale structure for 
each SN (note that it depends on i, or the redshift in 
general). The intrinsic scatter and the Poissonian large 
scale structure induced scatter add in quadrature, and 
both scale with 1 /Ni as expected i.e. they are both Pois¬ 
sonian in nature. The symbol quantifies the contri¬ 
bution to Cij from correlated/coherent large scale struc¬ 
ture fluctuations, and it does not scale inversely with the 
number of SNe. In principle, one should also add to Cij 
a term that describes systematic errors. We will discuss 
this briefly in but will leave it out for most of our 
discussion. We refer to Cij simply as the correlation ma¬ 
trix, to be distinguished from Cij which we call the mag¬ 
nitude covariance matrix i.e. Cij is the non-Poissonian, 
or coherent, part of Cij. 

The Poissonian large scale structure induced variance 
is simple to write down: 


/ ^Poiss. 



\ ^ r 

In 10 




( 6 ) 


which follows from eq. O- Here, Sd^iz) = SdL/dL{z) is 
the fractional fluctuation in luminosity distance due to 
departure from homogeneity (i.e. large scale structure) 
for an observed redshift z. Note that this is the fractional 
fluctuation at one point in the sky i.e. for one SN. The 
ensemble average ( ) is over realizations of the universe. 
The integration in the first line of eq. © is over the red¬ 
shift bin centered at Zi with width Azi. For a sufficiently 
narrow Azi, since ([^di, (-z)]^) typically varies slowly with 
z, the second line is a good approximation. 

The correlated/coherent component concerns Sd^ in 
different parts of the sky: 


C/j = 


\ ^ T 

In 10 


dzd^9dz'd^0' 

AziAiAzjAj 


{ddA^,e)SdAz',e')) ( 7 ) 


where we have added an extra argument 0 or O' to 
5dL to remind ourselves that the fluctuation depends on 
the redshift as well as the angular position. The an¬ 
gular integration is done over the area of the survey. 


and in general, the total area of the survey could be 
different at different redshifts, so we allow the area A 
to carry an index i or j. The redshift integration is 
done over the respective redshift bins. In cases where 
{Sdi^{z,6)SdL{z',0')) does not vary rapidly over the re¬ 
spective redshift bins (such as in the lensing contribu¬ 
tions to (5di), the expression above is well approximated 
by [5/lnl0p(AiAj)"^ / d^Od/^O'{5dL{zi,9)5dL{zj,0')). 

It is crucial that the Poissonian variance scales as I/W 
whereas Cij does not (eq. |Sj; this is justified in Appendix 
Naively one would expect Cij to be quite small es¬ 
pecially if the area Ai or Aj is large. However, for future 
surveys where Ni can be quite large, it is entirely possi¬ 
ble for the Poissonian term to be smaller than Cij , as we 
will see. 

It is also important to emphasize that Cij , even when 
i = j, is non-Poissonian, in the sense that i = j only tells 
us that we are looking at the same redshift bin, but 9 (or 
z) can still of course differ from 9' (or z'). In other words, 
when we talk about correlated/coherent/non-Poissonian 
fluctuations, we mean fluctuations that are correlated be¬ 
tween different SNe, at different redshifts or different an¬ 
gular positions or both. 

As far as we know, the importance of Cij has been over¬ 
looked in existing error forecasts for SN surveys. More¬ 
over, the only large scale structure contribution to the 
magnitude covariance matrix Cij (eq. |H1) that has been 
considered in the context of SN surveys is that from lens¬ 
ing, through the Poissonian term (crf”'®®-)^ (but see |23 
for preprints that partially overlap with ours concerning 
correlated lensing fluctuations). 


IV. THE LUMINOSITY DISTANCE 
FLUCTUATION AND ITS CONTRIBUTION TO 
THE MAGNITUDE COVARIANCE MATRIX 


Our goal in this section is two-fold. First, we write 
down an expression for , the fractional luminosity dis¬ 
tance fluctuation, that is accurate to hrst order in per¬ 
turbations. Second, we work out its two-point correlation 
and second moment as they show up in the magnitude 
covariance matrix (eq. E, [Hj and 0). The main results 
are eq. (tTHl) - CT . Readers not interested in details can 
skip ahead to m 


A. The Luminosity Distance Fluctuation 

There exist a large literature on the luminosity dis¬ 
tance in a weakly perturbed Friedmann-Robertson- 
Walker universe Sasaki ^3 derived the general ex¬ 
pressions using the optical scalar equations, and gave an 
explicit integration in the case of an Einstein-de-Sitter 
universe. More recently, Pyne and Birkinshaw [HO, 
using a different technique, derived explicit expressions 
for more general cases, including a non-flat universe as 
well as a universe with dark energy. In Appendix [nj we 
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present a derivation closely following that of min , but 
somewhat simplified. The end result differs slightly from 
13 , but only in terms that are subdominant compared to 
the gravitational leasing and peculiar motion terms, and 
so this difference has no impact on the main conclusions 
of this paper. The origin of the difference is explained 
in Appendix 10 Relegating the technical details to that 
Appendix, let us focus here on a somewhat heuristic but 
more intuitive derivation of the dominant terms, those 
arising from peculiar motion and leasing. 

Recall the following standard relation between the ob¬ 
served flux F and the intrinsic luminosity A, valid for a 
universe with or without inhomogeneities |3f)| : 


F{z) 


L 5VIq 
47r(l -I- zY 8Ae 


( 8 ) 


where is the proper area of the emitter, is the 
solid angle at the observer subtended by light rays from 
the emitter, and z is the observed redshift. The angular 
diameter and luminosity distances are respectively de¬ 
fined to be 


dA = \/5Ae/5^o , dL=dA{l + zY (9) 


We emphasize again that these expressions are equally 
valid in a homogeneous or an inhomogeneous universe. 

In a homogeneous universe, the above equations take 
the form 


F(Z): 

dA{z) 


L 

47r(l -f zY[dA{z)Y 
■ Xe/Y+z) 


Xe = x{z) 


dz'/H{z') 


^0 

Jl(z) = Ja(z)(1 -k zf 


( 10 ) 


where the bar on top of quantities reminds us that 
they are defined in a homogeneous universe. Here x is 
the usual comoving distance in a homogeneous universe, 
which is an integral over redshift of the inverse Hubble 
parameter F[, and Xe is the comoving distance to the 
emitter. 

Suppose we perturb such a universe by introducing pe¬ 
culiar motion, excluding for the moment other possible 
sources of fluctuations. This has two effects. First, z is 
Doppler shifted to a new value z: 


1-I-z = (1-I-z)(l-I-Ve • n - Vo • n) (11) 


OJ 

o 

a 

cd 

m 

• I—I 

t:3 

• (“H 

Vi 

O 

a 

• (“H 

s 



FIG. 1: A schematic Hubble diagram to illustrate the effects 
of peculiar motion. The solid line represents the luminosity 
distance in a homogeneous universe. Consider a SN repre¬ 
sented by the open square in such a universe. Suppose this 
SN is given a peculiar velocity moving away from the observer. 
This SN would then be displaced to the position of the solid 
square - note that both the redshift and the luminosity dis¬ 
tance are changed. The luminosity distance difference that 
we are interested in is 5dL{z) = dL{z) — dL{z), where z is the 
actual observed redshift, dL{z) is the actual observed luminos¬ 
ity distance, and dL{z) is the luminosity distance given by the 
solid line, at the same redshift z. The fractional luminosity 
distance ffuctnation that we study is = SdL{z)/dL{z). 


Peculiar motion causes (5Ho —> (5Ho(l — 2vo • n), which 
can be derived by performing boosts. 

Making use of eq. 0, CH) and m, we therefore have 


dL{z) = dL{z){l + 2ve-n-Vo-n) (13) 

It is important to keep in mind that the luminosity 
distance fluctuation that concerns us is (z) = [dL{z) — 
dL{z)]/dL{z), where and dL are at the same redshift 
z. Fig. n] illustrates this point. 

A Taylor series expansion of di tells us 


where Vg and vq are the peculiar velocities of the emitter 
and the observer respectively, and n is the unit vector 
from observer to emitter. This is accurate to first order 
in peculiar velocities. Note that the speed of light is set 
to unity. 

The second effect of introducing peculiar motion is to 
modify the angular diameter distance. Peculiar motion 
modifies SXIq but not 5Af. (to first order) in eq. 0 , 
leading to 

( 12 ) 


dL(z) = dL(z) ( 1 -b [1 -b -^][ve • n - Vo • n] ) (14) 

V <Xe J 

Note that the scale factor Oe, its derivative with respect 
to conformal time and the comoving distance Xe can be 
evaluated at either z or z: the changes in the expression 
above are of second order in the perturbations. 

Eq. lO and da together imply |37| 

(- 2 , n) = Ve • n- ^(ve-n-vo-n) (15) 

a'eXe 


dA{z) = dA{z){l -b Vo • n) 
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This is accurate to first order in peculiar velocities, ig¬ 
noring all other possible fluctuations. We have added 
the argument n to to emphasize the fact that 6di^ 
depends on direction (or angle) in addition to redshift. 

Gravitational leasing modifies the observed flux of an 
object without changing its redshift. The gravitational 
leasing magnification flensed^^uniensed jg worked out in 
many places e.g. [^1^ : 


magnification = 1-I-2 f dx-^ —(16) 
Jo Xe 

where (j) is the gravitational potential fluctuation, and 
is the Laplacian in comoving space. 

The gravitational leasing contribution to Sd^ is there¬ 
fore 


(17) 

Jo Xe 

In summary, the total peculiar motion and leasing con¬ 
tributions to dd[^ are 


Poiss., vel. \ 2 


r 5 1^ 

In 10 


X 






dJk k1 

(27r)3 k'^ 


P{k, a 


1 ) 


where P{k,a) and P{k±,a) represent the mass power 
spectrum at scale factor a and at wavenumber k and 
k = k± respectively, kz and k± denote the line-of-sight 
and transverse components of the wave vector (k± is the 
norm of the 2D vector kj^), D is the linear growth fac¬ 
tor and D' is its derivative with respect to conformal 
time, and any quantity with the subscript i is evaluated 
z = Zi- The symbol (erfstands for the Pois- 
sonian lensing term that is often studied: it gives the 
variance in convergence, up to a factor of [5/ In 10]^. The 
term ^(jP°iss. ,vei .^2 gjygg the variance in luminosity dis¬ 
tance due to peculiar motion, up to the same factor. Note 
that there is no cross velocity-lensing term - it vanishes 
because the lensing projection forces kz = Q. 

Likewise, from eq. o and eq. m, using the plane- 
parallel approximation (more discussion below), we And 


SdLiz,n) = Ve • n - — 

Xe La'Je 

{Xe-x)x^2 


(ve • n - Vo • n)(18) 


dx- 


Xe 


:V>(x) 


To reiterate: Ve and Vq are the peculiar velocities of the 
emitter and observer, and n is the line of sight unit vector 
pointing away from the observer (n here plays the role 
of 9 in eq. [7]); the comoving distance to emitter Xe, 
the scale factor at emission Oe and its derivative with 
respect to conformal time Og are evaluated at redshift z. 
One can see from above that for small Xe or at a low 
redshift, the peculiar motion term proportional to 1/xe 
becomes important, while at a large redshift, the lensing 
term (second line) is more important. A more rigorous 
derivation of Sd^ ^ together with an explanation of why 
other first order contributions can be ignored, is given in 
Appendix [Q 


B. From Sd^ to the Magnitude Covariance Matrix 


Our next task is to compute the second moment of 
for eq. and the two-point correlation of Sd^ for eq. 
0. This is carried out in detail in Appendix IDI Let us 
summarize the results here. 

From eq. © and eq. it can be shown that 


' Poiss.\2 / Poiss., lens\2 

) 


I / Poiss., vel. \2 
+ (f^i ) 


(19) 


[cr: 


Poiss., lens \2 _ 




5 

2 


In 10 


2 


( 20 ) 


r 

lo a 


2 


{x^ - x)x ' 

Xi 




;P{k±,a) 


_ ^lens I y^vel. 

ij ~ ^ij ^ij 


( 21 ) 


^iens _ 

^ij — 


r 5 

2 rmivi{xi,Xj) 

'{x^ - x)x 

InlO 2 
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where C*®"® and G™*' are the lensing and peculiar motion 
contributions to the correlation matrix G,,-. The window 
function for the lensing term IT/™® depends on the sur¬ 
vey geometry: the integration of 0 and 9' is over the 
survey area, which is allowed to depend on the redshift 
bin i or j for the sake of generality. The expression given 
for IT,*™® in terms of the Bessel function Ji assumes the 

‘-J 

survey spans a circular patch on the sky, with angular ra¬ 
dius or (i.e. Ai = 7r(6»““•)^ Aj = Tr{9f^^-y). 
The window function for the peculiar motion term IT™*' 
depends upon an integration over the comoving volumes 
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of the redshift bins: Vi and Vj. The expression given for 
in terms of the Bessel function and trigonometric 
functions assumes also that the survey spans a circular 
patch in the sky (see Appendix lUI for expressions appro¬ 
priate for more general geometries). Here and Ay^ 
are the widths of the respective redshift bins in comoving 
distance. Note also that IT™'' strictly speaking should 
have an imaginary part, which is odd in kz and so inte¬ 
grates to zero. For more complicated survey geometries, 
Wd™® and VF™'- can always be worked out from the def¬ 
initions above. 

The expressions in eq. m assume small angles or the 
plane parallel approximation. For the lensing term, this 
is acceptable since, as we will see, lensing is generally 
important only at high redshifts (z > 1) where practi¬ 
cal surveys of the future cover a sufficiently small area. 
However, for the velocity term, which is generally impor¬ 
tant at low redshifts where an ongoing survey such as the 
SNfactory covers a large area (almost half of the sky), it 
is useful to have an expression that does not assume the 
plane parallel approximation: 
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where n and n' are unit vectors pointing in the direction 
X and x', and k is the same for k. Since P{k, a = l)/k'^ 
does not depend on direction, one can replace Wj™' (k) 
by its average over the solid angle of k: 


w:;^-ik)=Y,i2£ + i) 

f^i d9'sine'Pe{cos9') 

/o 1 - cos6»“^’'' 
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dx 
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(23) 


d9 sm9Pi{ cos9) 
1 - cos9Y''^^- 


where ji is the spherical Bessel function, and j'^ is its 
derivative with respect to the argument, and Pi is the 
Legendre polynomial. Note that the survey area Ai = 
27r(l - cos6»““') (- 7r(6»““')2 only for small 6»““'). 

Eq. ll^ assumes a contiguous and circularly sym¬ 
metric survey geometry. Expressions for a more general 
geometry are given in Appendix I dI lea. ID 121 & Eni). 

Throughout this paper, all results presented come from 
using the exact expression for the velocity term (eq. l2‘Al ) 
rather than the small angle approximation for C™'' (eq. 
m)- It is an interesting question how well the small 
angle expression for the velocity term pproximates the 
exact expression. The answer is: surprising well even 


for surveys that cover large portions of the sky. This is 
discussed in the context of a concrete example in iYH 

Note that if 0™“- or equals tt (i.e. an all sky 

survey), only the monopole i = 0 term in eq. (I23II would 
be non-zero. More generally, if the SN survey covers 
large portions of the sky, only the low multipoles would 
be significant. 

Eq. (^^3^ , 1^^^^ , and are the main results 

of this section. Their derivation from eq. m is given in 
Appendix lUl They are to be substituted into eq. 0 to 
obtain the magnitude covariance matrix. 

One subtlety: eq. (PI - ^ assume the vq term in 
eq. m can be removed, because the observer’s pecu¬ 
liar motion is fairly well-known from the dipole of the 
microwave background |^. However, when the SN is 
sufficiently close to the observer, care must be taken to 
take into account the correlation between Vg and Vq in 
eq. m- This is discussed in Appendix [HJ For realistic 
SN surveys, we have checked that the correlation makes 
a negligible difference. An alternative error estimate is 
to allow Vq to be a random variable and include its con¬ 
tribution to the variance of the magnitude fluctuations - 
doing so would increase the errorbar. 

It is worth emphasizing that the above expressions are 
derived using linear perturbation theory. For instance, 
the velocity terms in eq. (PI - (|22l assume linear per¬ 
turbation growth. This is acceptable for C™*' (eq. 
and [221) since the velocity integrals are dominated by 
large scale modes for relevant survey areas. This is 
more questionable for the Poissonian velocity contribu¬ 
tion 1201)) since this depends 

on the one-point (small scale) velocity dispersion. How¬ 
ever, as we will see in ESI the velocity contribution to 
the Poissonian term is relatively unimportant compared 
to the intrinsic contribution, and it is not necessary to 
get it exactly right. To be precise: we will use the linear 
power spectrum for the velocity terms (Poissonian and 
otherwise) throughout this paper. We use the transfer 
function from |4(l| . assuming a scalar spectral index of 
ns = 0.95, a Hubbble constant of h — 0.7, a matter den¬ 
sity of Dm = 0.27, a baryon density of fib = 0.046, and a 
normalization of as = 0.8 |^. (Our fiducial cosmology 
is always a flat cosmological constant dominated universe 
with Ode = 0.73; see the end of lUIl i 

For the lensing integrals in eq. (O - (EU , it is not 
uncommon to use the nonlinear mass power spectrum in 
the integrands assume the results are a good 

approximation even on small angular scales (a rigorous 
justification can be found in ^3)- Whether one uses the 
linear or nonlinear power spectrum does not matter much 
for C*®"®' (eq. [23 ) which is dominated by large scale 
modes anyway for realistic survey areas. For the Poisso¬ 
nian lensing term j-gq^ (201)) the 

nonlinear mass power spectrum would provide a more ac¬ 
curate estimate. Moreover, it is possible that the result is 
further enhanced (beyond that predicted by the standard 
nonlinear mass power spectrum) by the adiabatic con¬ 
traction of halos due to the cooling of baryons, presence 

























of MACHOs, etc. This will be further discussed in 
Throughout this paper, the lensing integrals are always 
done using the nonlinear power spectrum prescribed by 
[ 4 ^ (with suitable enhancement for the Poissonian term 
to account for plausible effects of the baryons; see below). 


V. PUTTING EVERYTHING TOGETHER - 
FROM LUMINOSITY DISTANCE 
FLUCTUATIONS TO ERROR FORECASTS 

Our goal in this section is to make error forecasts by 
putting together the results from the last 3 sections, en¬ 
capsulated in eq. lO for the Fisher matrix, eq. © for 
the magnitude covariance matrix , and eq. ira - ra 
for explicit expressions for components of the magnitude 
covariance matrix in terms of the power spectrum (with 
a small modihcation to eq. cni; see below). It will prove 
illuminating to hrst study separately the velocity and 
lensing contributions to the covariance matrix, and their 
effects on the parameter errors. This is done in ESI 
and IVRl In WU we make error forecasts for several 
ongoing/planned/proposed SN surveys, and discuss the 
relative importance of the velocity- and lensing-induced 
luminosity distance fluctuations. 


A. Peculiar Motion 

Let us consider a survey with only one redshift bin 
with a width of Az = 0.05. We are interested in the ve¬ 
locity contributions to the Poissonian and non-Poissonian 
terms in the magnitude covariance matrix (eq. [H]; with 
one redshift bin, the covariance matrix is just a num¬ 
ber Cii). Note that at the low redshifts we consider in 
this section, the lensing contributions are negligible. The 
coherent/non-Poissonian velocity term Cjf' (eq. ED) is 
plotted as solid lines in Fig. [3 It is shown as a function 
of the mean redshift of the survey. In other words, Az is 
fixed at 0.05, but the central redshift is allowed to vary, 
and Cif' is shown as a function of that central redshift. 
The various solid lines span from top to bottom a survey 
area of 1000 to 41000 square degrees ^3- ^ contigu¬ 
ous circuliarly symmetric geometry is assumed, except 
for the second solid line from the bottom: it consists of 
two circular patches one centered at the north and the 
other south, with a total area of 20000 square degrees, in 
other words, it has a galactic cut of about ±30°. 

For comparison, we show as a dashed line the Poisso¬ 
nian velocity term ™^ )^/iV (eq. and [201): 

with N = 300 SNe. Note that this term is independent 
of the survey area. We also show as dotted lines the in¬ 
trinsic Poissonian term (cr™*'' )^/iV, for cr™*''- = 0.1 (lower 
line) and cr™*''- = 0.15 (upper line). 

The above numbers are chosen for a reason: the Nearby 
Supernova Factory (SNfactory) is a survey of about 300 
SNe that covers roughly half of the sky (^ 20000 square 
degrees), and is centered around z ~ 0.055 spanning a 



FIG. 2: Various contributions to the magnitude covariance 
matrix Ci\ (eq. |^) for a low redshift SN survey characterized 
by one redshift bin with Az = 0.05. The five solid lines show 
the coherent/correlated velocity term Cif' (eq. 1211 1 as a 
function of the mean redshift z, for a survey area of, from top 
to bottom, 1000, 5000, 20000, 20000b and 41000 (full sky) 
square degrees. A contiguous circularly symmetric geometry 
is assumed, except for the case of 20000b, which has the 20000 
square degrees split into two patches, one centered at the 
north pole and the other south (i.e. a galactic cut of about 
±30°). Note that the SNfactory (SNf) has a mean 2 of 0.055 
(2 = 0.03 — 0.08), and a total area of about 20000 sq. degrees 
(denoted by the open square). The dashed line shows the 
Poissonian velocity contribution (eq. 1201 1 to the magnitude 
covariance matrix: '— 3 gg 

note that it is independent of the survey area, but depends on 
the mean 2 . The dotted lines show the contribution from the 
intrinsic magnitude scatter: ((t 1“‘'' )^/A, where N = 300, and 
equals 0.1 for the lower line and 0.15 for the upper line. 
Lensing contributions to Cii are negligible at these redshifts. 


width of Az ~ 0.05. (Other similar surveys include the 
Carnegie Supernova Project and LOTOSS, see 0.) To 
be precise, henceforth, whenever we discuss the SNfac¬ 
tory, we assume a geometry that coincides with the next 
to bottom solid line of Fig. H (The moral of the bot¬ 
tom three solid lines is that neither the precise area nor 
the precise geometry matters much, as long as the survey 
covers a significant fraction of the sky and does not have 
many holes or edges.) One can see that for these parame¬ 
ters, the coherent fluctuation term C™*' is larger than the 
Poissonian (intrinsic and velocity) terms, and becomes 
even more important if the survey were done at a lower 
redshift. The Poissonian velocity term (dashed line), on 
the other hand, is always subdominant compared to the 
intrinsic term |4q - this is probably the reason for the 
common perception that peculiar motion can be ignored 
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as part of the error budget for SN surveys. This percep¬ 
tion is incorrect because it ignores the coherent velocity 
fluctuations quantified by Cjf' (solid lines). 

It is useful to understand qualitatively why the dif¬ 
ferent contributions to Cn depicted in Fig. Intake the 
values they do. The contribution from the Poissonian 
intrinsic scatter is the simplest: 0.1^/300 or 0.15^/300 
giving 3.3 x 10“® or 7.5 x 10“^ (dotted lines). The con¬ 
tribution from the Poissonian velocity term is also easy 
to understand. At low redshifts, the term 
[201) is roughly 2.17 x v/(cz) where v is the typical pecu¬ 
liar velocity (~ 300 km/s) and cz is the Hubble flow. For 
instance, at z = 0.055, this amounts to ^ 0.04, 

and therefore ™^')^/300 ~ 5 x 10“® (dashed line). 

For the non-Poissonian velocity term (eq. 1221 1. 

let us focus on the case corresponding to the SNfactory, 
with a total area of 20000 square degrees and a mean 
redshift of z = 0.055 (the lowest black line). The large 
survey area means that the window function 
is dominated by the low order multipoles (eq. 1231 1. 
Let us consider the monopole £ = 0, which picks out 
k ~ 0.005 h/Mpc corresponding to a mean distance of 
X ~ Ay; ~ 200 Mpc/h. The integral over power spectrum 
(second line of eq. [221 ) can therefore be approximated 
by 47rfcP(fc)/(27r)^ evaluated at fc ^ 0.005 h/Mpc, giving 
roughly 4 (Mpc/h)^. The prefactors in the first line of eq. 
(P|l equal ~ (2.17)^ x (1/0.055)2 x (0.5/3000)2( h/Mpc)^, 
where we have made use of the fact that D' is roughly 
half the inverse Hubble radius ~ 0.5/(3000Mpc/h) (re¬ 
call that the speed of light is set to one). Putting all 
these together yields ~ 2 x 10“"^. 

A low redshift survey such as the SNfactory provides 
an important anchor for surveys at higher redshifts in 
that it helps determine the zero-point M (eq. PP). As 
we will see, combining high redshift SN surveys with a 
low redshift survey such as the SNfactory often reduces 
the error on the equation of state of dark energy by a 
factor of about 2. It is therefore important to ask: to 
what extent does peculiar motion, particularly coherent 
peculiar motion, increase the projected error on M from 
a survey like the SNfactory? 

Fig. [^provides the answer. The dotted lines show the 
errorbar on M (keeping all other parameters hxed |4^1 
from a survey of 300 SNe that spans z = 0.03 — 0.08, ig¬ 
noring peculiar motion i.e. only the intrinsic magnitude 
scatter is taken into account: the upper dotted line is for 
an intrinsic scatter of = 0.15, and the lower dot¬ 
ted line is for cr™*"'- =0.1. With only the intrinsic scatter 
taken into account, the error on AT is independent of sur¬ 
vey area. The solid lines show the same, except this time 
including peculiar motion induced fluctuations. As be¬ 
fore, the upper line of the pair uses cr'"*’'' = 0.15 and lower 
line uses cr'"*''- =0.1. (At these redshifts, other sources of 
large scale structure fluctuations such as lensing are neg¬ 
ligible.) For a survey like the SNfactory (^ 20000 square 
degrees), one can see that peculiar motion increases the 
error on AT by about a factor of 2, depending on the 
intrinsic scatter assumed. This result makes good sense 



FIG. 3: The zero-point (M in eq. ^) rms error as a func¬ 
tion of survey area (keeping all other parameters fixed). The 
survey redshift coverage is fixed: z = 0.03 — 0.08, and the 
number of SNe is 300. The upper pair of solid lines allow for 
the effects of peculiar motion, while the lower pair of dotted 
lines do not. Within each pair, the upper line uses an intrin¬ 
sic scatter of = 0.15 and the lower one uses cr'"'*''' = 0.1. 
Note that the SNfactory (SNf) covers half of the sky, which 
is about 20000 square degrees. 

because we can see from Fig. |21that including the coher¬ 
ent velocity contribution C/®*' raises the total magnitude 
covariance by a factor of 3 — 4. The lesson: peculiar ve¬ 
locity has a significant impact on the determination of 
the SN zero-point from a low redshift anchor. 

Note that in all our computations of the coher¬ 
ent/correlated velocity term we use the exact ex¬ 

pression that allows for large angles (eq. [211 )• We find 
that using the plane parallel approximation (eq. m) 
leads to an underestimate of by only about 10%, 
even for a survey with high sky coverage like the SNfac¬ 
tory. 

Fig. 121 might give one the impression that one is better 
off moving the low redshift anchor to a higher z where 
the peculiar motion induced magnitude fluctuations are 
smaller, largely because the ratio of peculiar velocity to 
Hubble flow is smaller. However, to measure the equation 
of state Wpivot accurately, it is advantageous to have a 
long lever arm in redshift. In other words, a large redshift 
span (from the low redshift anchor like the SNfactory to 
a high redshift SN survey) is preferable - recall that the 
discovery of cosmic acceleration comes from comparing 
the low redshift part (z < 0.1) of the Hubble diagram 
with the high redshift part (z ^ 1). From this point 
of view, it is not immediately obvious that moving the 
low redshift anchor to a higher z actually helps. Fig. 01 
addresses this question. 
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FIG. 4: Marginalized error on Wpivot as a function of the 
mean redshift of a low z survey of 300 SNe, with A 2 = 0.05. 
The upper panel assumes an intrinsic magnitude dispersion 
of = 0.15, and the lower panel uses = 0.1. Here, 
it is assumed the low 2 ; survey is combined with a high 2 
SNAP-like survey whose parameters are kept fixed (see Table 
EJ. The solid lines show the error on Wpivot when peculiar 
motion induced fluctuations (dominated by the coherent ones) 
are taken into account: the upper solid line is for a low 2 
survey area of 1000 square degrees and the lower one is for 
20000 square degrees. The dotted line shows the same when 
peculiar motion is ignored, hence the result is independent of 
survey area. For comparison, the horizontal dashed line shows 
the error on Wpivot if one uses only the data from the high 2 
SNAP-like survey. Gravitational leasing is taken into account 
in all cases above. A prior of = 0.03 (flat universe) is 
assumed. 


The solid lines show the marginalized error on Wpivot 
when a high z SNAP-like survey of 2000 SNe (see Table 
^for details) is combined with a low z survey of 300 SNe, 
width A 2 = 0.05 and a mean 2 as shown on the x-axis 
(with a geometry like that used for the SNfactory, see 
Fig. 121). A prior of rms = 0.03 is assumed. The 
two solid lines are for a low 2 survey area of 1000 (up¬ 
per) and 20000 (lower) square degrees respectively. They 
allow for all sources of magnitude fluctuations we have 
discussed: peculiar motion, gravitational lensing and in¬ 
trinsic scatter (the upper panel is for cr™*"'- = 0.15 and 
the lower panel uses cr™*"'- = 0.1). One can see that the 
optimal mean redshift for the low 2 survey is about 0.08. 
However, as long as the mean redshift is above 0.06 or 
so, the precise redshift does not appear to matter much. 
The upturn of the error at 2 < 0.06 is due to coherent 
peculiar motion as discussed before (Fig. [2J|. The flat¬ 
ness of the solid lines at 2 > 0.06 is due to the rough 
cancellation of two opposing effects mentioned above: a 


higher redshift (for the low redshift anchor) is good for 
suppressing peculiar motion induced fluctuations while a 
lower redshift is useful for creating a long lever arm. The 
SNfactory has a mean redshift which is sufficiently close 
to optimal that it is probably not worth moving it to a 
higher redshift where observations are more challenging. 

The dotted line shows the error on rcpivot for exactly 
the same set up as above but with peculiar motion ig¬ 
nored (the result is independent of the survey area of 
the low 2 survey). One can see that by ignoring peculiar 
motion, one might reach the erroneous conclusion that 
the SNfactory should be moved to a lower redshift. We 
emphasize that it is the coherent peculiar motion that 
matters here - the Poissonian velocity induced fluctua¬ 
tions are simply too small to be of consequence (see Fig. 
m. For comparison, we also show in Fig. with a hor¬ 
izontal dashed line the error on Wpivot if one uses only 
the data from the high 2 survey. The difference between 
the dashed line and the solid/dotted lines illustrates the 
benefit of having a low 2 anchor. 

Rather similar conclusions are reached about the op¬ 
timal redshift of the low 2 anchor when we examine the 
impact on the error of Wa, or when we combine the low 
2 survey with other high 2 surveys e.g. those in Tabled 


B. Gravitational Lensing 

Here, we are interested in the relative importance of 
the Poissonian and non-Poissonian lensing induced fluc¬ 
tuations, The latter has been 

largely overlooked in error forecasts for SN surveys. As 
we will see, unlike the case of peculiar motion, it is the 
Poissonian term that is important for lensing. Because 
this term is sensitive to small scale fluctuations, as men¬ 
tioned in 3TVT^ we always use the nonlinear mass power 
spectrum to compute all lensing quantities. 

To facilitate comparison with previous treatments of 
the Poissonian lensing fluctuations, we introduce one 
modification to our previous expressions: replace eq. ra 
by 


Wi ) = 2 ( 0 -^ 


Poiss., lens\2 


r + K' 


(24) 


The boost factor of 2 appears to be necessary to roughly 
reproduce the results of 0 who used a halo approach 
instead of followin g a power spectrum approach like we 
do. (The results of [ij] are about 2 — 3 times hi ghe r than 
ours without the boost factor.) In other words, [l^ effec¬ 
tively used a different power spectrum from ours. Phys¬ 
ically, such a boost could arise from non-gravitational 

t sics (recall that the nonlinear power spectrum we use 
arises purely from gravitational instability). For in¬ 
stance, the adiabatic contraction of dark matter halos 
due to the radiative cooling of baryons can enhance the 
power spectrum on small scales |48| . There could even 
be a signihcant population of MACHOs which can lens 
the SNe (see e.g. 0 )- On the other hand, there is 
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FIG. 5: Various lensing contributions to the magnitude co- 
variance matrix (eq. [^). Upper panel: the solid lines show 
the non-Poissonian lensing term (7*™° (eq. 1211 1 as a function 
of z = Zi for a survey of 1, 5 and 15 square degrees (top to 
bottom); the dotted line shows the Poissonian lensing term 
2 X (eq. & 1201 1. for N = 100; for reference, 

the dashed lines show the expected contribution from intrin¬ 
sic scatter for = 0.15 (upper line) and 0.1 

(lower line). Lower panel: the normalized lensing correla¬ 
tion across redshifts i.e. as a function of 

z = Zi, for three different values of Zji 0.045, 0.95 and 1.95. 

at least some evidence suggesting that Cold Dark Mat¬ 
ter models might overpredict the amount of small scale 
structure. The precise boost factor is therefore a bit un¬ 
certain. The choice of 2 is somewhat arbitrary, but it 
seems prudent to include some enhancement of power 
due to well-motivated physical effects ^3. Ultimately, 
high redshift SNe themselves will tell us what the right 
level of Poissonian lensing fluctuation is. In any case, the 
boost leads to a more conservative errorbar. Throughout 
this paper, eq. (El is used in place of eq. m when 
making error forecasts. 

It is interesting to note that if the small scale power 
spectrum has a shape similar to that seen in N-body 
simulations the integral for |^ ) 

is dominated by fc ~ 10 h/Mpc. This coincides with 
the scale where the adiabatic contraction of halos due 
to baryon cooling, in some sense the most plausible non- 
gravitational effect, is expected to become important 
On larger scales, gravity is almost certainly the only sig¬ 
nificant shaping force of large scale structure. This is 
why no boost factor is necessary for the non-Poissonian 
lensing term C]™® (eq. ETl ). 

Fig. 0 shows the relative importance of the Poisso¬ 
nian and non-Poissonian lensing terms. The upper panel 
shows with solid lines the non-Poissonian lensing term 


FIG. 6: Similar to Fig. except that the survey area is 24 
square degrees, and N = 700, numbers that are motivated by 
JEDI (think of N as roughly the number of SNe per Az of 
0 . 1 ). 


where we have divided the redshift into bins of 
Az = 0.1 each The term C]®"®' depends on sur¬ 

vey area, and we show from top to bottom the result for 
1, 5, and 15 square degrees. The dotted line shows the 
Poissonian lensing term 2(crf°‘®'*'’where one can 
think of fV = 100 as the number of SNe in a redshift bin. 
For comparison, we show with dashed horizontal lines 
the Poissonian term due to intrinsic scatter (cr'"*’' )^/100, 
with the upper line using cr™*"'- = 0.15 and the lower line 
using =0.1. The number of N = 100 is somewhat 
arbitrary but it roughly corresponds to the number of 
SNe per redshift bin of 0.1 in a survey like SNAP, which 
also has an area of around 15 square degrees (Table P). 

What can we conclude from Fig. O’ From the upper 
panel, one can see that lensing begins to be of compara¬ 
ble importance to the intrinsic scatter when the redshift 
climbs above 1 or so. Moreover, unless the survey area 
is quite small (less than ^ 1 square degree), the non- 
Poissonian lensing term is small compared to the Poisso¬ 
nian lensing term at redshifts where they matter. From 
the lower panel, we can see that there is a signihcant 
amount of lensing induced correlations between different 
redshifts. However, because C)®"® is small to begin with 
(compared to /N), we do not expect these cross- 

redshift correlations to hugely impact dark energy errors 
from SN surveys. We have verified this to be the case. 

The above conclusion about the relative importance of 
non-Poissonian and Poissonian lensing terms is subject to 
changes in survey parameters, however. Fig. Elis similar 
to Fig. El except that the number of SNe per redshift bin 
is increased to 700, as motivated by the ambitious SN sur- 
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vey JEDI. (The area is also slightly increased to 24 square 
degrees from the 15 square degrees of SNAP, though that 
has a minor impact relatively speaking.) One can see that 
the non-Poissonian and the (boosted) Poissonian leasing 
terms are here of comparable importance at the redshifts 
where they matter {z > 1). 


C. Error Forecasts: Peculiar Motion versus Lensing 


Survey 

No. of SNe 

Area (sq. deg.) 

Redshift Dist. 

DES 

1900 

40 

0.2 - 0.8 fiat 

ESSENCE 

200 

12 

0.2 - 0.8 flat 

JEDI 

14000 

24 

1 

o 

SDSSII 

200 

250 

0.05 - 0.35 flat 

SNAP 

2000 

15 

O 

i 

SNfactory 

300 

20000 

0.03 - 0.08 flat 

SNLS 

600 

4 

0.2 - 0.8 flat 


TABLE I: A summary of survey parameters we study in this 
paper. Each set of parameters are supposed to mimic, but 
not necessarily exactly match, actual SN surveys that go by 
these names Q. For instance, the (flat) redshift distributions 
almost certainly differ from the actual ones - they are chosen 
for simplicity. For JEDI and SNAP, the redshift distribu¬ 
tion is taken from M, scaled to the appropriate number of 
SNe. For some surveys not included in this table, such as the 
Carnegie Supernova Project and CfA Supernova Program at 
low redshifts and the LSST and Pan-STARRS at high red- 
shifts, see Q. 

We are finally ready to make error forecasts for sev¬ 
eral examples that approximate ongoing/proposed SN 
surveys. They are described in Tabled Note that the 
adopted parameters by no means exactly match those 
of actual surveys that go by those acronyms, but they 
should be close. Also, except for the SNfactory (see Fig. 
13, we assume the area listed for each survey is over a 
contiguous region, and, for simplicity, consists of a circu¬ 
lar patch on the sky. (We will have more to say later on 
the implications of the exact survey geometry.) See 0 
for further details on these surveys, and others we have 
not worked out explicitly. Our procedure, simply put, 
is: we put eq. (1^ and eq. ® - 123 for the large 
scale structure induced magnitude fluctuations into eq. 

and eq. for the Fisher matrix, and obtain the 
relevant dark energy errors. 

For readers interested in the bottom line, Fig. [7| & 
121 are in some sense the most important figures of this 
paper. They show the degradation in the marginalized 
error on Wpivot for several high redshift SN surveys, each 
used in conjunction with the SNfactory as a low redshift 
anchor. Degradation refers to the fractional increase in 
error. It is defined as JrCpivot— 1, where 
is the rms error when only the intrinsic scatter is in¬ 
cluded, and SwpWot is the rms error when all sources 
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FIG. 7: The degradation in the marginalized error on the 
equation of state Wpivot for several high redshift SN surveys, 
each used in conjunction with the SNfactory as a low redshift 
anchor. (The error on iCpivot is marginalized over Wa, Dde 
and M; because of our choice of the scale factor pivot, this 
is identical with the error on a constant equation of state, 
marginalized over fide and M; see am) Degradation refers to 
the fractional increase in the rms marginalized error on Wpivot 
due to large scale structure induced fluctuations i.e. degrada¬ 
tion = Jn)pivot/(5uip“vot “ li where dtCpTot allows for only the 
intrinsic scatter while Juipivot takes into account both the in¬ 
trinsic scatter and large scale structure. The division of each 
histogram into black and white regions show how much of the 
degradation is due to lensing (black) and how much due to pe¬ 
culiar motion (white). The upper panel uses 0.1 for the intrin¬ 
sic scatter while the lower panel uses 0.15. In all cases, 

except JEDI, the lensing degradation is completely dominated 
by Poissonian lensing fluctuations. For JEDI, the black-filled 
portion is the degradation due to Poissonian lensing fluctua¬ 
tions and the black-hatched portion is the additional degra¬ 
dation due to non-Poissonian/coherent lensing fluctuations. 
A prior of Jflde = 0.03 (fiat universe) is assumed. 


of fluctuations (intrinsic -I- lensing -|- velocity) are in¬ 
cluded (full histogram, including both black and white 
portions), or when only intrinsic + lensing fluctuations 
are included (black portion of the histogram). In other 
words, the white (black) portion of the histogram tells 
us the degradation due to peculiar motion (lensing). In 
all cases other than JFDI, the lensing fluctuations are 
completely dominated by the Poissonian ones (the ones 
that are customarily considered). For JFDI, the black¬ 
filled portion represents the degradation when intrinsic 
-I- Poissonian-lensing fluctuations are included, and the 
black-hatched portion represents the additional degrada¬ 
tion when non-Poissonian/coherent lensing fluctuations 
are taken into account. The error on Wpivot here is ob¬ 
tained by marginalizing over Wa, fide and M (as discussed 
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FIG. 8: Same as Fig. Qexcept that no prior on Ude is assumed, 
though the universe is still flat. 

in m and the Appendix HI this is identical to the error 
on a constant equation of state, marginalized over fide 
and M). A prior of rms (5flde = 0.03 is assumed for Fig. 
[3 while no such prior is used for Fig. jS] A flat universe 
is assumed in both cases. The upper panel uses an in¬ 
trinsic scatter of size cr™*"'- = 0 . 1 , and the lower panel 
uses = 0.15. An important conclusion: the general 
trend is for surveys that have a higher statistical power 
to suffer more degradation, and in most cases, peculiar 
motion dominates the degradation (white portions of the 
histograms). 

Fig. El shows the error contours for the example of 
a SNAP-SNfactory combination. The dotted contours 
include the intrinsic scatter only, while the solid con¬ 
tours include both the intrinsic scatter and large scale 
structure (velocity -I- leasing) fluctuations. The dashed 
contours are for the case where some level of systematic 
error = 0.02(1 -|- z)f2.7 [^) is included on top of 
intrinsic plus large scale structure fluctuations. One can 
see that a systematic error of this size degrades the de¬ 
termination of rcpivot to a degree that is comparable to 
large scale structure (see the end of this section for more 
discussion). Note that this figure uses cr'"*''- =0.1; for 
^mtr. _ 0 25 ^ constraints would be worse. 

For readers interested in further details, Table UH dis¬ 
plays error forecasts for all interesting parameters for 
various combinations of surveys, assumptions and priors. 
There are several interesting points to be made. 

• The prior of = 0.03, which is often assumed 
in the literature, makes quite a difference to the 
projected errors on the equation of state Wpivot and 
its slope Wa (eq. El), especially for the smaller 


FIG. 9: Projected (68.3%) errors for SNAP -|- SNfactory 
for three different combinations of sources of noise: intrin¬ 
sic scatter only (dotted), intrinsic scatter -|- large scale struc¬ 
ture induced fluctuations (solid), and intrinsic scatter -|- large 
scale structure -|- systematic error (dashed). The upper panel 
shows the errors for tide and Wpivot marginalizing over the 
derivative Wa = —dwfda and the zero point M, assuming a 
flat universe but no prior on tide- The lower panel shows the 
errors for Wa and Wpivot marginalizing over the dark energy 
density tide and the zero-point M, assuming a flat universe 
and a prior of (Aide = 0.03. Here, w = Wpivot + Wa(apivot — a), 
where ttpivot is chosen to make the errors on uipivot and Wa 
uncorrelated (see 23- Note that Swa ~ 25w' , where w' is an¬ 
other common parametrization of the evolution of the equa¬ 
tion of state: w' = dw/dz. In both panels, (j‘“*''' = 0.1 is 
used. The degradation in errors due to large scale structure 
is primarily due to coherent peculiar motion and Poissonian 
lensing fluctuations (see Fig. EI^ISl- The systematic error 
is assumed to be = 0.02(1 -|- z)/2.7 IboII . It appears 

the degradation due to large scale structure and that due to 
systematic error are comparable. 

surveys - those with a narrower redshift range and 
fewer SNe. 

• The addition of a low redshift survey like the SNfac¬ 
tory generally improves the errorbars significantly, 
the more so for surveys with less statistical power. 

• As mentioned in m the error on Wpivot, marginal¬ 
ized over the slope Wa, dark energy density fide and 
the zero-point M, is exactly equal to the error on 
a constant equation of state, marginalized over fide 
and M (this is the result of an optimally chosen 
scale factor pivot). This is relevant especially for 
the smaller surveys - they produce only weak con¬ 
straints on the variation of w, which is why the 
error on w generally projected for these surveys as¬ 
sumes no variation with z. Our marginalized error 
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on Wpivot can be compared directly against these 
forecasts. Note, however, that our error on fide, 
marginalized over Wpivot, Wa and M, is generally 
larger than the error on fide, marginalized over a 
constant equation of state and M. 

• Several surveys have a rather similar redshift cover¬ 
age of 0 - 0.2-0.8: DES, ESSENCE and SNLS. On 
their own, without the addition of the SNfactory, 
one can see that the large scale structure fluctua¬ 
tions make little difference to the errors (compare 
’air with ’intr’). This is because peculiar motion is 
only important for z <0.1 and lensing is important 
only for 2 : > 1. It is with the addition of the SNfac¬ 
tory (which is useful because it improves the errors 
quite a bit) that peculiar motion has an impact. 

• Large scale structure fluctuations generally have a 
larger impact (fractionally) on the errors for larger 
surveys, a point already illustrated in Fig. [7|and|Hl 

• For the most part, we have so far focused on sta¬ 

tistical errors. This is in part because the level 
of systematic errors is uncertain, depending on the 
precise capabilities and operational details of the 
respective experiments. To get some idea of how 
systematic errors might change our conclusions, we 
try two different prescriptions for the systematic 
error: one from [T^ = 0.02 x {z/1.7), and the 

other from = 0.02 x {1 + zj^.l) (see 

on redshift binning). They behave similarly at high 
0 , but the latter implies a larger systematic error 
at low z. We apply these to SNAP -I- SNfactory 
in Table QH with the former prescription labeled 
as ’sys’ and the latter labeled as ’sys2’. The gen¬ 
eral conclusion is that systematic error contributes 
to a degradation that is comparable to large scale 
structure. This is illustrated in Fig. compare 
the increase in error for Wpivot from the dotted to 
solid contours against that from the solid to dashed 
contours. 


VI. DISCUSSION 

Let us summarize the main lessons. 

• Large scale structure induced magnitude fluctua¬ 
tions have a significant impact on dark energy mea¬ 
surements from a whole array of ongoing or future 
SN surveys. For instance, the degradation, due to 
large scale structure, in the error for the equation 
of state Wpivot ranges from 10% to 60% depending 
on surveys and assumptions (Fig. [7||SI). It appears 
difficult to measure from SNe alone the equation of 
state to better than about 7 — 10% (depending on 
assumptions), unless one has a survey considerably 
more ambitious than SNAP -|- SNfactory (see Fig. 

13 and Table mil. 


• Of all possible large scale structure fluctuations, the 
dominant ones are due to peculiar motion and grav¬ 
itational lensing. Peculiar motion is important at 
2 < O.I (through a low redshift anchor such as the 
SNfactory) while lensing dominates at z >1. The 
impact of peculiar motion is mainly through co¬ 
herent/correlated large scale flows (Fig. [SJ while 
the impact of lensing is mainly through Poissonian 
fluctuations (Fig. [3). The Poissonian fluctuations 
can be reduced by increasing the number of SNe, 
while the coherent ones can only be suppressed by 
increasing the survey area. When a high redshift 
(z >0.1) survey is combined with a low redshift an¬ 
chor (z < O.I), we find that peculiar motion mostly 
dominates over lensing as a source of error. 

• What does the above mean for survey designs? As 
has been emphasized in the literature (e.g. 0 ), a 
low redshift anchor such as the SNfactory is very 
useful for reducing the eventual dark energy errors 
from a high redshift SN survey. For such a low 
redshift survey, one might have hoped to reduce 
the coherent peculiar motion induced fluctuations 
by either increasing the survey area or moving it 
to a higher redshift. All else being equal, neither 
will improve appreciably the precision on dark en¬ 
ergy determination. The SNfactory already covers 
half of the sky; going to full sky will not reduce 
the errors significantly. For instance, combining 
SNAP with an all-sky version of SNfactory instead 
of the half-sky one that we have been assuming, 
the marginalized error for SwpWot would improve 
by only about 1%. Furthermore, moving SNfac¬ 
tory to a higher redshift, while useful in reducing 
peculiar motion induced fluctuations, shortens the 
lever arm that the combination of a high redshift 
survey and a low redshift anchor offers. The net 
effect is that moving the low redshift anchor to a 
higher redshift actually does not reduce the error 
on the equation of state rcpivot appreciably (Fig. 0)). 

• How about survey designs for a high redshift sur¬ 
vey? For a high redshift survey (considered on 
its own) that does not extend beyond 2 : ~ I 
(e.g. DES, ESSENCE and SNLS all cover roughly 
2 ^ 0.2 — 0.8), neither peculiar motion nor lens¬ 
ing constitutes significant sources of errors. The 
only way to reduce dark energy errors is to in¬ 
crease the number of SNe, and suppress systematic 
errors. The precise survey area is of little impor¬ 
tance for such a survey, as long as it is not too 
small (too small meaning I square degree or less, 
see Fig. 0. (Of course, peculiar motion does play a 
role in the eventual errors once one combines such 
a high redshift survey with a low redshift anchor, 
which as emphasized above, is generally a good 
idea.) For a high redshift survey that extends be¬ 
yond 2 ~ I, gravitational lensing becomes a non- 
negligible source of errors. But because lensing’s 
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Survey 

Jwpivot (prior) 5wa (prior) 

6M (prior) 

Jwpivot (no pr.) 

3wa (no pr.^ 

(5Dde (no pr.) SM (no pr.) 

DES -b SNf (all) 

0.085/0.079 

0.69/0.57 

0.021/0.018 

0.17/0.14 

4.9/3.5 

0.47/0.33 

0.027/0.023 

DES -b SNf (intr) 

0.073/0.061 

0.60/0.46 

0.013/0.009 

0.13/0.09 

4.5/3.0 

0.44/0.29 

0.018/0.012 

DES (all) 

0.11/0.01 

1.4/1.0 

0.064/0.044 

0.36/0.25 

17/12 

1.4/0.9 

0.23/0.15 

DES (intr) 

0.11/0.10 

1.4/1.0 

0.063/0.042 

0.35/0.23 

17/11 

1.3/0.9 

0.22/0.15 

ESSENCE-bSNf (all) 

0.094/0.087 

1.5/1.1 

0.024/0.021 

0.36/0.26 

14.0/9.6 

1.4/1.0 

0.043/0.033 

ESSENCE-bSNf (intr) 

0.084/0.075 

1.4/1.0 

0.017/0.011 

0.34/0.23 

13.6/9.1 

1.3/0.9 

0.038/0.025 

ESSENCE (all) 

0.19/0.15 

4.2/2.9 

0.20/0.13 

1.11/0.76 

54/37 

4.2/2.9 

0.69/0.47 

ESSENCE (intr) 

0.18/0.14 

4.1/2.8 

0.19/0.13 

1.1/0.72 

52/35 

4.1/2.7 

0.68/0.45 

JEDI -b SNf (all) 

0.041/0.032 

0.46/0.41 

0.013/0.011 

0.042/0.033 

0.76/0.60 

0.042/0.033 

0.017/0.013 

JEDI -b SNf (intr) 

0.032/0.021 

0.41/0.32 

0.010/0.007 

0.032/0.022 

0.60/0.40 

0.033/0.022 

0.011/0.008 

JEDI (all) 

0.046/0.035 

0.52/0.46 

0.018/0.014 

0.048/0.036 

0.99/0.73 

0.052/0.039 

0.025/0.018 

JEDI (intr) 

0.041/0.028 

0.50/0.42 

0.017/0.013 

0.042/0.028 

0.90/0.60 

0.046/0.031 

0.024/0.016 

SDSSII + SNf (all) 

0.12/0.10 

5.7/4.0 

0.039/0.031 

1.6/1.1 

101/69 

15/10 

0.092/0.067 

SDSSII -b SNf (intr) 

0.105/0.081 

5.4/3.6 

0.031/0.021 

1.5/1.0 

99/66 

15/10 

0.085/0.057 

SDSSII (all) 

0.17/0.12 

8.0/5.4 

0.081/0.056 

2.3/1.6 

145/97 

21/14 

0.20/0.14 

SDSSII (intr) 

0.17/0.12 

7.8/5.2 

0.078/0.052 

2.2/1.5 

143/95 

20/14 

0.20/0.13 

SNAP + SNf (all) 

0.069/0.058 

0.57/0.52 

0.018/0.015 

0.079/0.063 

1.5/1.2 

0.09/0.07 

0.022/0.019 

SNAP + SNf (intr) 

0.054/0.038 

0.53/0.46 

0.011/0.008 

0.06/0.04 

1.32/0.88 

0.078/0.052 

0.014/0.009 

SNAP + SNf (all-bsys) 

0.075/0.067 

0.60/0.55 

0.018/0.016 

0.090/0.077 

1.8/1.5 

0.11/0.09 

0.022/0.020 

SNAP -b SNf (intr-bsys) 

0.061/0.049 

0.56/0.50 

0.012/0.008 

0.073/0.056 

1.6/1.2 

0.10/0.08 

0.015/0.010 

SNAP + SNf (all-bsys2) 

0.081/0.074 

0.63/0.59 

0.020/0.018 

0.101/0.089 

2.0/1.7 

0.12/0.10 

0.025/0.023 

SNAP + SNf (intr-bsys2) 

0.071/0.062 

0.59/0.55 

0.015/0.013 

0.086/0.073 

1.8/1.5 

0.111/0.093 

0.019/0.016 

SNAP (all) 

0.101/0.078 

0.70/0.62 

0.038/0.028 

0.121/0.087 

2.6/1.9 

0.14/0.10 

0.066/0.046 

SNAP (intr) 

0.096/0.069 

0.68/0.60 

0.037/0.026 

0.110/0.074 

2.4/1.6 

0.122/0.082 

0.063/0.042 

SNLS -b SNf (all) 

0.089/0.084 

0.98/0.75 

0.022/0.019 

0.24/0.18 

8.2/5.7 

0.81/0.56 

0.032/0.026 

SNLS -b SNf (intr) 

0.077/0.069 

0.89/0.64 

0.014/0.009 

0.21/0.14 

7.9/5.3 

0.78/0.52 

0.025/0.017 

SNLS (all) 

0.14/0.12 

2.5/1.7 

0.114/0.078 

0.64/0.44 

31/21 

2.4/1.7 

0.40/0.27 

SNLS (intr) 

0.13/0.11 

2.4/1.6 

0.11/0.07 

0.62/0.42 

30/20 

2.4/1.6 

0.39/0.26 


TABLE II: Marginalized la errors for different combinations of surveys and assumptions. The survey parameters are taken 
from Table E](SNf here stands for the SNfactory). The descriptions in parentheses, ’ah’, ’intr’, ’sys’, ’sys2’ refer to the sources 
of errors that are included in the forecasts - ’all’ means including all sources of random fluctuations intrinsic-scatter -|- lensing 
-I- velocity, ’intr’ means including only the intrinsic scatter, and ’sys’ and ’sys2’ mean systematic error (see text for details). 
The errorbar on each parameter is obtained by marginalizing over all the other parameters (altogether, there are four: Wpivot, 
Wa, fide and the zero-point M; see anj. Those errorbars with the description ’prior’ assume a prior of rms dflde = 0.03, while 
those denoted with ’no pr.’ assumes no such prior. For each entry, we give two numbers in the form x/y, where x assumes the 
intrinsic scatter has a size of = 0.15 and y assumes = 0.1. Note that in cases where the errorbars are sufficiently 
large, the Fisher matrix analysis likely breaks down. The above numbers are rounded-off at two significant figures. A flat 
universe is assumed throughout. 


impact is mainly through the Poissonian fluctua¬ 
tions it introduces, increasing survey area (such as 
for SNAP) is not really necessary. The only in¬ 
stance in which a case can be made for increasing 
survey area is JEDI, which has a sufficiently small 
Poissonian error (due to its large number of SNe) 
that coherent/correlated lensing fluctuations actu¬ 
ally play a role (see Fig. [7|and|Hl). 

Our investigations in this paper naturally raise a num¬ 
ber of questions and issues, some of which we address 
briefly here, and some require further research. 


1. Perhaps the most natural and interesting question 
is whether current constraints on dark energy, which typ¬ 
ically come from some combination of high and low red- 
shift SNe {z > 0.1 and < 0.1), are already affected by 
peculiar motion. The short answer is: not very much. 
This is because the current number of low redshift SNe 
used (typically several lO’s) is sufficiently small that the 
Poissonian error (due to simply intrinsic scatter) is quite 
a bit larger than the coherent velocity error. This can be 
inferred from Fig. |21 raising the dotted lines by a factor 
of ^ 10 (due to dropping the number of SNe from 300 
as in the figure to ~ 30) means the Poissonian intrinsic 
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scatter constitutes a larger source of error than coherent 
peculiar motion (solid lines). Note that this argument 
assumes the existing low redshift SNe are selected from 
a large area of the sky e.g. (so that it is the lowest 
few solid lines of Fig. 13 that is relevant). The conclusion 
could be quite different if this assumption does not hold. 
We urge SN experiments to clearly state the survey areas 
of their different samples (especially the low z samples) 
when publishing their results. 

2. In our forecasts for a selection of representative SN 

surveys (Table we have assumed simple geometries - 
a contiguous circular region on the sky for the high z 
surveys, and two separate patches (one in the north and 
one in the south) for the SNfactory. Realistic surveys are 
bound to be more complicated in shape, possibly with 
many holes or gaps. The exact geometry affects the size 
of the coherent fluctuations but not the Poissonian ones. 
Since for the most part the only coherent fluctuations 
we need worry about are those due to peculiar motion 
which is important only at low z’s, it is mainly the exact 
geometry of something like the SNfactory that concerns 
us. It is therefore worth repeating our calculations for 
the actual geometry of the SNfactory, including possible 
extra gaps for instance. For this purpose, we have given 
sufficiently general expressions for the relevant window 
functio ns (W /r" and ) in eq. (EIJ and (HU, and 
in eq. and in Appendix lUl Note that gaps 

almost always increase the importance of coherent fluctu¬ 
ations because of the introduction of high k modes. For 
the high redshift surveys that extend beyond z ~ 1, it 
would be useful to check that a realistic survey geome¬ 
try does not make the coherent lensing fluctuations much 
more important (though we do not expect this to happen, 
as long as the survey area exceeds ^ 1 square degree). 
For high z surveys that stay within the redshift range 
O.I <z <1, neither velocity nor lensing fluctuations are 
expected to be important (unless the number of SNe is 
much larger than what has been considered), and so the 
survey geometry has a relatively minor impact. 

3. We have largely ignored internal motion in our dis¬ 
cussions of velocity induced fluctuations. By internal mo¬ 
tion we mean the motion of the SNe within galaxies, for 
instance due to the virialized motion of the SN progen¬ 
itors, or even due to the orbital motion of the SN itself 
within the binary system that is its progenitor. Such 
motion could contribute to the overall peculiar velocity 
of the SNe. Ignoring internal motion is partially justi¬ 
fied by the fact that in practice the redshifts are assigned 
based on the redshifts of the host galaxies. (We thank 
the referee for emphasizing this to us.) However, internal 
motion could still in principle modify the apparent lumi¬ 
nosity. The important point to keep in mind is that such 
internal motion is not expected to be correlated between 
SNe in different galaxies, and so our calculation of the 
coherent/correlated velocity fluctuations remains valid. 
Internal motion can certainly increase the Poissonian ve¬ 
locity fluctuations. However, typical virialized motion is 
of the order of a few hundred km/s, similar to the typical 


large scale flow velocity, and so the Poissonian velocity 
fluctuations remain subdominant (Fig. [SJ. Also, the or¬ 
bital motion internal to the binary progenitor is too slow 
to be of significance (see e.g. jl^). 

4. Our main focus in this paper is on statistical errors: 
from intrinsic scatter and from large scale structure in¬ 
duced fluctuations. We have investigated the effect of 
systematic error in some simple examples (see Table IH) 
entries for ’SNAP -h SNf (all -f sys)’ and ’SNAP -H SNf 
(all -I- sys2)’, and the associated discussion at the end 
of WU see also Fig. 0 - we show that even in the 
presence of systematic error of the assumed magnitudes, 
large scale structure fluctuations remain a non-negligible 
source of errors for dark energy measurements. It would 
obviously be useful to investigate this further and explore 
a wider range of systematic errors suitable for each SN 
experiment. 

5. An implicit assumption in our calculations is that 

the redshift measurements of low z SNe are sufficiently 
accurate for us to worry about their peculiar motion in 
the first place. Existing low z measurements typically 
report an accuracy of 6z ~ 0.001 — 0.002 (e.g. [^1. 

The spectral instrument of the SNfactory has a resolu¬ 
tion of 1200, corresponding to Sz ~ 0.001. (Note that 
the actual redshift accuracy is likely to be better than 
the instrumental spectral resolution, so this is a conser¬ 
vative estimate.) Translating into velocities, we are talk¬ 
ing about a velocity of 300 — 600 km/s, or a magnitude 
fluctuation of 6m ~ 0.04 — 0.08 (for z ~ 0.05). This is 
still smaller than the intrinsic magnitude scatter that we 
assume: cr™*"'' = 0.1 or 0.15. The redshift uncertainty 
adds to the Poissonian scatter, and our range of 0.1 to 
0.15 can be thought of as accounting for this possibil¬ 
ity already. It is important, however, that the redshift 
measurements do not suffer from a systematic bias (that 
affects all SNe in the same way). From Fig. |3 it can be 
seen that a systematic bias of 5z ~ 0.0003 or 100 km/s 
(at z ~ 0.05) would have a comparable effect as coherent 
peculiar motion. 

6 . An interesting question is: to what extent can 
corrections be made for the velocity and lensing fluctu¬ 
ations? For instance, one could imagine using galaxy 
weak lensing maps to correct for the magnihcation of 
SNe. This has been shown to be not viable, or not suf¬ 
ficiently accurate to be useful, by [s^- This is because 
galaxy weak lensing maps typically tell us the magnihca¬ 
tion on scales larger than are relevant for the Poissonian 
part of SN lensing. (These maps can be useful for cor¬ 
recting the non-Poissonian/coherent part of SN lensing, 
but this part of lensing is not very important for most 
SN surveys anyway.) More recently it was argued by 
that corrections for (the Poissonian part of) SN lensing 
can be made by modeling foreground galaxies as isother¬ 
mal spheres or generalizations thereof. One should keep 
in mind that the Poissonian lensing fluctuations are sen¬ 
sitive to structures on relatively small scales (fc >10 
h/Mpc), and a smooth halo proHle does not necessarily 
capture all the relevant fluctuations. For instance, in the 
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case of strong lensing, substructures are often invoked to 
explain the observed flux ratios |^. Most of the high 
z SNe will not be strongly lensed, but a similar lesson 
applies here. Moreover, using the foreground galaxies to 
make a magnification correction inevitably involves as¬ 
sumptions about bias: how galaxies trace mass. As can 
be seen from Fig. El /'Jn ~ 0.01 for SNe at z ~ 1.5 
and N = 100 (per Az of 0.1, as appropriate for SNAP 
for instance). For the lensing correction to be useful, it 
should therefore satisfy two criteria: first, the correction 
should be more accurate than 0.1 (in magnitude) per SN; 
second, it should not introduce a systematic bias that is 
larger than 0.01 (in magnitude; magnitude fluctuation 
^ magnification fluctuation). Even if the first can be 
achieved, the second seems challenging. 

How about corrections for velocity fluctuations? Here, 
the situation is slightly different: what needs to be cor¬ 
rected is the coherent part (i.e. large scale), not the Pois- 
sonian part, of the fluctuations. Roughly speaking, we 
need to know the low order multipoles of the peculiar flow 
at the redshift of e.g. the SNfactory (z ~ 0.055). One 
option is to use peculiar velocity surveys (such as from 
the SNe themselves), but it should be kept in mind that 
to disentangle the Hubble flow from say the monopole, 
one needs a survey that has the same sky coverage as 
the SNfactory, but is deeper. This requires considerable 
resources. Another option is to use the galaxy spatial 
distribution as a guide, i.e using mass conservation to 
relate peculiar velocity to the galaxy overdensity. Such 
a procedure of course suffers from the uncertain biasing 
relation between galaxies and mass. Note also that one 
needs a galaxy survey that is deeper than the SNfactory 
to define the correct mean galaxy density. (For a recent 
paper that examines the peculiar motions predicted by 
the PSCz survey, see it focuses on peculiar flows at 
slightly lower redshifts than we need.) Whether either 
option allows us to take out the effect of bulk flows to 
sufficient accuracy is a question we would like to address 
in the future. 

7. Another natural and interesting question is: to 
what extent can the noise here, due to lensing and pe¬ 
culiar motion, be viewed as a useful signal? In the case 
of lensii^ the issue is discussed in several recent papers 
In general, it is difficult for SNe to be 
competitive with galaxies as the sources for weak grav¬ 
itational lensing experiments. Consider for instance the 
measurement of the convergence power spectrum: the 
shot-noise in the case of SNe is (cr™*"'')^/(4n) (the fac¬ 
tor of 4 comes from 5m ~ 2k where k is convergence), 
while the shot-noise in the case of galaxies is {a'^Y/{2n) 
(the factor of 2 comes from the use of two components of 
shear to estimate k). Here the intrinsic magnitude scat¬ 
ter cr™*’'- is roughly 0.1 — 0.15, the shape noise cr^ is about 
0.3, and the surface density n is approximately 0.04 per 
square arcminute for SNe (taking numbers from SNAP), 
while n ~ 30 per square arcminute for a typical weak 
lensing galaxy survey. The shot-noise from SNe is sim¬ 
ply too big compared to that from galaxies. Nonetheless, 


the lensing of SNe is free from certain systematic errors 
that might affect the lensing of galaxies, such as intrinsic 
alignment, and so the SN method still provides a useful, 
though not terribly stringent, consistency test. 

How about the SN peculiar motion as a signal? For 
some of the earlier work on this issue, see e.g. 

1^ l63| . Among the different methods for measuring pe¬ 
culiar velocities (see for a review), SNe la constitute 
the most accurate distance indicator on an object by ob¬ 
ject basis. For instance, SNe la yield distances with an 
error of ^ 5 — 7%, while Tully-Fisher distances are typi¬ 
cally uncertain at the 15 — 20% level. On the other hand, 
Tully-Fisher galaxy catalogs (e.g. [b^) typically have 

significantly more objects than SN surveys, and there¬ 
fore have perhaps more statistical power. Yet, SN sur¬ 
veys might suffer less from systematic errors that seem 
to have plagued at least some Tully-Fisher galaxy cata¬ 
logs (e.g. | 6 J), and SN surveys generally go deeper. It 
remains an interesting question to what extent compet¬ 
itive cosmological constraints can be obtained from the 
peculiar motion of SNe. We hope to explore this in the 
future. 
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APPENDIX A: ON ERRORBARS 

Let us briefly describe how to go from the Fisher ma¬ 
trix to errorbars by giving some examples. Labeling the 
parameters Wpivot, Wa^ f^de and M as pi to p 4 , the (rms) 
errorbar on Wpivot marginalized over everything else with 
no prior is given by If a prior on Ode is desired, 

say with rms Sflde = 0.03, the marginalized errorbar on 
Wpivot is equal to ^/[{F + G)“i]ii, where G is a diago¬ 
nal matrix with all zero-entries except for G 33 = 1/0.03^. 
On the other hand, the error on Wpivot keeping everything 
else fixed is y/l/Fu, which is generally smaller than the 
marginalized error -^[F~ ']ii- For instance, in this pa¬ 
per, we have occasion to study yO/i/jd, the error on M 
keeping everything else fixed. 

It is also useful to elaborate the nature of the pivot 
in the parametrization: w{a) = Wpivot + iCa(apivot — a)- 
Suppose we start with choosing Opivot = Ij and from 
the Fisher matrix, after appropriate marginalization over 
Ode and M (with or without prior), arrive at a co- 
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variance matrix for the errors on Wpivot and Wa of the 
form: Qn = Q12 = Q21 = (Swpi^otSwa) and 

Q22 = (<5^0) ■ For instance, if we marginalize over both 
fide and M with no prior, one can obtain Q from invert¬ 
ing the full Fisher matrix: Qn = Q12 = 

and Q22 = [F~^]22- In general, the choice of Opivot = 1 
would not lead to a vanishing Qi2. Choosing a different 
Opivot from the original Opivot = 1 corresponds to a simple 
linear transformation on the parameters Wpivot and Wa- 
Our goal is to choose a new Opivot such that the errors 
on uipivot and Wa are uncorrelated. This can be accom¬ 
plished by choosing Opivot = 1 + Qi2lQ22- With this 
choice, the resulting error on rcpivot (marginalized over 
Wa) becomes Qu — Q12IQ22, while the error on Wa 
(marginalized over Wpivot) remains the same i.e. \JQ22- 
It can be shown that this marginalized error on tCpivot is 
exactly the same as the error on w if dwjda were fixed 
(to be zero, say) [sj- This is easiest to see if there are 
only two parameters involved: Wpivot and Wa- With our 
choice of Opivot, Q is turned diagonal and so is its inverse 
(the Fisher matrix). In such a case, marginalization over 
Wa is equivalent to fixing Wa (to be zero for instance). It 
is straightforward to generalize this to the case of more 
parameters (e.g. including Ode and M). In other words, 
the error on Wpivot, marginalized over Wa, Ode and M, 
is exactly the same as the error on a constant equation 
of state, marginalized over Ode and M. This is a useful 
fact to know because some experiments quote errors on 
w assuming a vanishing dw/da. 

For completeness, the following relation is useful: the 
background dark energy density scales with redshift as 
p (X (1 -I- exp [— 3 wo 2;/(1 -|- z)] for the 

parametrization w{a) = iCpivot + ^0(1 — a). This follows 
from energy momentum conservation d{pa^) = —Pd{a^), 
and P = w{a)p. 


APPENDIX B: THE MAGNITUDE 
COVARIANCE MATRIX I - GENERALITIES 

Our goal is to derive eq. © for the magnitude co- 
variance matrix Cij, as well as eq. 0 and 0 for its 
individual components. 

The average magnitude rrii from a redshift bin i is: 

= (Bl) 

^ a 

where W is the number of SNe in that bin, denotes 
the magnitude of individual SN labeled by a (instead of 
the averaged magnitude as in rm), Qai equals 1 if SN a 
falls within bin i and vanishes otherwise. The summation 
of a is over all SNe in one’s survey. Note W = ®ai- 
From above, we have 

Cij = {SmiSmj) = ^((5m°(5TOg)0„0&j (B2) 

* ^ a,b 


The summation over a and b can be split into two terms: 
a = b and a ^ b. The a = b term gives: 

-^^^^(((5m°)2)0„,0„j = ^^(((5m°)2)0„i (B3) 

^ ^ a * a 


The magnitude variance of an individual SN, (((5m°)^), 
has two contributions, one from intrinsic variations and 
the other from large scale structure induced fluctuations. 
We usually treat the intrinsic part as independent of a. 
Therefore, we have 


1 ,_. ('_intr.'i2 i / Poiss.'i2 

^ (B4) 


W; 


where (cr™*"' )^ is intrinsic, and (crf°‘®® )^ is due to struc¬ 
tures. The latter is defined to be 


1 




(B5) 


which can be approximated by taking the continuum 
limit: 



(B6) 


This last expression makes use of the fact that Sm^ = 
(5/ln 10)<5di loJ' IF® large scale structure fluctuations (eq. 
[P), and that ([^di,]^) depends on redshift but not angular 
position. 

Eq. IIB4II reproduces exactly the Poissonian terms in 
eq. 0 , while eq. 03 matches eq. 0 . 

It remains to study the a ^ b contribution to eq. ifR^ . 
Once again, making use of i5m° = (5/lnl0)i5d^ for the 
large scale structure induced fluctuations, we can approx¬ 
imate this by taking the continuum limit: 


1 


'^{Sm°Jml)QaiQbj 

a^b 


(B7) 


r 5 1 ^ 

In 10 


dzd^Odz'd^d' 

AziAiAzjAj 


{5a,{z,e)5d,{z',e')) 


There is no need to zero out the z = z' and 9 ~ 9' con¬ 
tributions in the above integrals because they are van¬ 
ishingly small. The expression above matches the non- 
Poissonian term Cij in eq. 0 and O- 

If one wishes, one can add a term, in addition to those 
in eq. 03 and 03, to {SrriiSmj) accounting for sys¬ 
tematic errors. 


APPENDIX C: LUMINOSITY DISTANCE IN A 
PERTURBED FRW UNIVERSE 

Our aim is to derive an expression for the luminosity 
distance fluctuation that is accurate to first order, and 
show that among all first order terms, those in eq. 
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dominate in realistic applications. The derivation follows 
closely the one given by Pyne & Birkinshaw mui The 
final expression agrees with that of [13 in the case of 
an Einstein de-Sitter universe. It also agrees with the 
earlier result from Sasaki who also gave an explicit 
expression for an Einstein de-Sitter universe. It disagrees 
slightly with |17| for a more general universe, but only 
in terms that are subdominant compared to the lensing 
and velocity terms kept in eq. The origin of this 

difference is discussed below. 

For simplicity, we assume a flat universe. The metric 
in a convenient gauge, the conformal Newtonian gauge, 
takes the form 

+ 2(/))d7?2 + {l- 20) (Cl) 

i 

= a"[-(l + 20)d7;2 + (1 - 20)(dx^ + x^dn^)] 

where 0 is the scalar perturbation which is also the grav¬ 
itational potential, ry is the conformal time, 0 ( 77 ) is the 
scale factor, x® is the Cartesian comoving coordinate 
{i = 1,2,3), X is the radial comoving distance, and dfl^ 
is the usual angular part of the metric. We ignore vec¬ 
tor and tensor fluctuations and assume zero anisotropic 
stress. 

We will take the path of first working out the angular 
diameter distance, and then using the general relation 
eq. to obtain the luminosity distance. 

The null geodesics can be most easily worked out by 
ignoring the factor of in ds^, since it is just an overall 
conformal factor. For such a metric {ds^/a?), the affine 
connection components are 

r%o = <^,., = r%, = r'oo = 0,z(C2) 

r Oj — ^ 3 ^ — ^jk4^^i 

We use Latin indices to denote the spatial components. 
Note that even though we work out the null geodesics in 
the rescaled metric ds^/a^, all subsequent manipulations 
(i.e. subsequent to eq. 1(151 1 are done with the original 
metric ds^. 

Let A be the affine parameter, and let us split the pho¬ 
ton path into background and first order pieces: 

77 (A) = 77 (A) -I- 5ri[X) , x*(A) = x*(A) -I- 5x\X) (C3) 

The zero-order background solution is just a straight 
line 


77 (A) = rio — X, x*(A) = Xrd (C4) 

where 770 is the conformal time today, and rd is a unit vec¬ 
tor pointing away from the observer. Note the slightly 
strange choice of the sign for A: it increases as one 
traverses backward in time along the path of the pho¬ 
ton. The choice is perfectly acceptable since the geodesic 
equation is invariant under a sign flip for A. 


The first order geodesic equation can be integrated to 
give 

f'X nX 

Sr]{X) = 2 dX'cj) + 2 dA'(A - A')0.r, (C5) 

Jo Jo 

i5x*(A) = —2 f dX'{X — X')[di — rdn^dj]<j) 

Jo 

-2 [ dA'(A-A')0,„n* 

JO 


where we have chosen the boundary conditions (5x*(0) = 
577 ( 0 ) = d5xVdA(0) = 0, with A = 0 denoting the loca¬ 
tion of the observer. The condition on the photon di¬ 
rection at the observer {dx^/dX{0) = dx’'/dX{0) = rd) 
implies d6rj/dX — 20 vanishes at the observer by nullness 
of the geodesic. Note for objects with Latin (spatial) 
indices, we are being cavalier about their placement up¬ 
stairs or downstairs: e.g. rii = rd is the /i = 7-th compo¬ 
nent of n^. 

Note that the photon 4-momentum in the original met¬ 
ric ds^ is given by 


k. = l^ 
dX 


(C6) 


The factor oil/a3 is necessary because the affine param¬ 
eter A was chosen in the rescaled metric ds^ jc?. See [^ 
for a proof. 

We are interested in a cone of light emanating from 
the observer back towards the emitter i.e. the observer is 
located at the tip of the cone. Suppose the center of this 
cone points in the (fiducial) direction n (at the observer). 
A photon moving along this direction follows a geodesic 
described by the above equations and (EH), with 
n = n. Consider a photon moving along the surface of 
this cone i.e. its path direction at the observer differs 
slightly from the fiducial direction, say n = n-|-5n. Such 
a photon follows a geodesic also described by eq. E3 
and ESI), but with n = n -I- 5n. 

Our strategy is to work out x^(A) for these two paths 
which tells us how the light cone gets deformed as A 
moves away from zero (the observer) and allows us to 
compute the angular diameter distance. 

Let’s define the deviation w^\ 


w^(X, n) = x^ 
_ dx^(X,n) 
~ 077^ 


(A, n -I- 5n) 
6n^ 

n—n 


— x^(A, n) 


(C7) 


where we have added the argument n to x^ to remind 
ourselves that these are geodesics that start out at the 
observer with a particular direction. 

We are interested in the area spanned by at 
the emitter. More precisely, it is the area in a two- 
dimensional space perpendicular to both (for the fidu¬ 
cial ray) at the emitter and the four velocity of the emit¬ 
ter. Let us define a projection operator M that 
projects onto this space: 


w 






(C8) 
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Suppose for simplicity that n = (0,0,1) i.e. points in 
the direction. Then (5n points in the — x^ plane. 
One can parametrize it by 

i5n = (e sintj}, e cosiIj, 0) (C9) 


where ip can take any value from 0 to 27r, and e basically 
defines the size of the light cone at the observer’s position 
(though the observer sees a light cone with a different size 
if she/he has a non-zero peculiar velocity; see below). 

Now, consider how varies as we vary ijj: 


< = 




9a;'" (A, n) 


9ni 

9a;'"(A, n) 
dv? 


esin^/ (CIO) 

I 

ecos'tjj 


In the two-dimensional space where lives, this traces 
out an ellipse. In other words, as a function of ip 

can be rewritten in the form sin^ {ip — tpo) + cos^ {ip — 
ipo) where p and q are the lengths of the major and minor 
axes. The area of such an ellipse is irpq. Relating p and 
q back to the quantities in eq. (IHTOll . it can be shown 
that the area of the light beam at the emitter is 

SAe = (Cll) 

\ 1/2 

-(K]iKm]2)") 


where we have defined 

(C12) 

n=n 


n=n 

Eq. icriTi and icTl . together with the geodesics 
worked out earlier, allow us to compute the light beam 
area at the emitter. One fact simplifies the calcula¬ 
tion a bit: the quantities [w^ 2 {wx^i ]2 and 

{w^i[w±fP \2 can be computed accurately to first or¬ 
der by just using the zero-order H^^,. The zero-order 
for our choice of coordinates is simple, it is di¬ 
agonal: = diag(0,1,1,0). To see this, note that 

the only possible place where we might care about the 
first order part of H^i, is when it is contracted with 
the zero order part of [9x'"(A, n)/9n^] or [dx''{\, n)/9n^], 
which is simple: (0,A, 0, 0) or (0, 0,A, 0) (see eq. 1(141 1. 
This means the only part of the first order piece of 
that we care about is H^i or H^ 2 - Moreover, 
in the quantities of interest such as = 

^[dx’'{X, n) /dn^]H°'pldx^{X, n) /dn^]gfj_a, suppose we 
are interested in the contribution from the first order 
piece of the first - then only the zero-order pieces 
of the other terms need to be considered, and to zero- 
order, both and are diagonal (and = 

diag(0,1,1,0)). The upshot is that the only first order 
piece of that we need to worry about is H‘^ 2 , 




FIJ2 


9x'"(A,n) 

9a;'" (A, n) 
dv? 


H ^2 or It can be shown explicitly that the first or¬ 
der corrections to these vanish. Further discussions can 
be found in [T3 |. 

Putting everything together, we can substitute the 
zero-order and the geodesic solution from eq. 
and HU into eq. iim . and get: 


[w]_]i=X + 2[ dA'(A-A')</.3 (C13) 

Jo 

pX pX 

-2 / dX\X- X')X'(l)^ii-2 / dV(A-V)0,r? 

Jo Jo 

[w‘]_]i = [w ]_]2 = —2 f dX'{X — X')X'(j)^i 2 
Jo 

[ujj_]2 = A + 2 / dA^(A — X')(j)^^ 

Jo 

pX pX 

-2/ dA'(A-A')A>,22-2 / dX'{X- X')(p,r, 

Jo Jo 


One useful tip in deriving the above is to remember that (p 
in the integrands is not only an explicit function of A", but 
also an implicit function of n (i.e. d(p{X')/dn‘^ = X'(pp{X') 
to first order). The result above can be put into eq. icTfiTi 
to obtain the area at emission: 


5Ae=TTe^al{l-2(Pe)Xl (C14) 



where Ae is the affine parameter at the emitter, and all 
quantities with a subscript e are evaluated at emission. 

To find the angular diameter distance, we also need 
to know the solid angle of the beam according to the 
observer. If the observer does not have peculiar motion, 
the answer is simple: from eq. (I(19ll . the solid angle is 
obviously tte^. If the observer has peculiar motion, the 
solid angle can be obtained by performing a simple boost: 

Jflo = 7re^(l - 2?;i^) (C15) 

where Vq is the line of sight component of the observer’s 
peculiar velocity i.e. dx^/dg. 

With the above two equations, all ingredients are in 
place to write down the angular diameter distance and 
therefore also the luminosity distance (eq. |H1): 

dL{Xe) = aeXe{l + Zef (^1 - (pe + (C16) 

-2 J ^ [-(p^3 + (p^r,] 

We also need an expression for 1 + Ze, the observed red- 
shift of the emitter. This can be obtained from the ratio 
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where is the photon momentum in 
eq. Kihll and is the 4-velocity of the emitter/observer: 
= dx^/dr = [(1 — (/)/a, v/a]. Therefore 


1 + z(Xe) 


1 


!-[</>- - 2 / ' ^,r,d\' 

Jo 


(C17) 


where [Q]q for some quantity Q is defined to be Qe — Qo, 
the difference at emission and observation. 

Ultimately, we are interested in the luminosity distance 
fluctuation Sdi^ = [^^(Ae + JAe) — (iL(Ae)]/(iz,(Ae), where 
dL is the luminosity distance in an unperturbed universe, 
and (5Ae is chosen such that 1 + z(Ae + SXe) = 1 + z{Xe) 
with z being the redshift in the unperturbed universe (see 
Fig. nj. 

From eq. iITTitIi . it can be shown that to first order 


1 + z{Xe + (5Ae) — [1 + z{Xe 


( 


1 + 


a 

a 

.3ie 


{SXe - 2 / (/>dA' - 2 


^^,rjdX') 



(CIS) 

A')(/.„dA') 


To derive this, it is useful to remember that the scale 
factor for affine parameter Ae+(5Ae is equal to a evaluated 
at the time fj{Xe + (5Ae) + Sr](Xe), accurate to first order. 
The scale factor evaluated at f]{Xe) is of course equal to 
l/[l + z(Ae)]. 

Similarly, from eq. 1C16II . one can show that 


dbiXe + <5Ae) — a[77(Ae)][l + z(Ae + <5Ae)]^Ae (C19) 




((5Ae - 2 / </dA' - 2 / (Ae - X')(j),r,dX') 


a. 


— (j)e + 'Cq “I-^-2 


dX 


,Ae-A' 


9,3 i- 9,Vi 




Choosing 6Xe such that l+z{Xf,+5Xe) = l+z(Ae) in eq. 
(IClSIl . and noting that (/^(Ae) = a[? 7 (Ae)][l + z(Ae)]^Ae, 
we finally have from eq. lITTmll : 


5d^ {z) =vl-(l)o-2cj)e-2 dX'cj)^^ (C20) 


AgUg 


[</ - v% + 2 / dX'cj)^ 


/o 


dX' 


— dX --- [ 9,11 + 0,22] 

/■Ae \ 

+ 2 / dW \ (hr. 


where z is the observed redshift in the perturbed uni¬ 
verse. All quantities with subscript e are to be evaluated 
at Ae; evaluating them at Ae + SXe instead would only 
make a difference to to second order. 


It is useful to rewrite the expression above in the fol¬ 
lowing form: 


SdL ( 2 , n) = Ve • n-— (ve • n - Vo • n) (C21) 

Xe La'-le 


Jo 


Xe 


1 


a 


Xe La'Je 


(</e ~ /’o) + 2 


dx 


Hx) 


10 


0 Xe 

’‘■,X2izxirix)-2rdx^4.'(x) 

Xe Jo Xe 


Xe ia'i 


a 


dx(t>'{x) 


which is accurate to first order in fluctuations, and we 
have switched from n pointing in the direction to al¬ 
low it to point anywhere (and we have relabeled n —> n). 
Note that </' = </>,,,, (j)" = (p^rjv ^^.d a' = The argu¬ 
ment X (radial comoving distance) of (j){x) is supposed 
to remind us that </> in general depends on position, as 
well as, implicitly, time i.e. under the line of sight in¬ 
tegrals, (j){x) should be evaluated at position x and at 
a time when a photon is supposed to reach x- Sym¬ 
bols with subscript e denote evaluation at the time of 
photon emission and symbols with subscript 0 denote 
evaluation today. One manipulation we have done is to 
replace + (j ),22 in the integrands by ~ 4>,33 and 
use d/dx^ = d/dX' X-d/dr]. Also, we have replaced Ae by 
Xe where Xe is the zero-order comoving distance to the 
emitter. Doing so is justified because fluctuations in Xe 
make a difference to the final expression above only to 
second order. To be precise: 


Xe= f dz'/H{z') (C22) 

Jo 

where z is the observed redshift in the perturbed universe 
(same z as the argument of Sd^^ in eq. Eni). 

The above expressions (eq. IC^ & Eni) are consis¬ 
tent with eq. (5.16) of Sasaki who gave an explicit 
expression for an Einstein-de-Sitter universe (adopting 
Sijo = 0 in ^3)- Pyne & Birkinshaw’s da eq. (32) can 
be most readily compared with our eq. (1(12011 . The two 
are identical, except that 0 did not have the last term 
on the right hand side of our expression: 


2 



dX 


,Ae-A' 


(C23) 


It takes a little bit of work to trace back the origin of 
the difference. Everything up to (and including) the ex¬ 
pressions for dL and the redshift, eg. and (irTTTIi . 

appears to be consistent with (e.g. our eq. 1(1141 
& Em are consistent with their eq. 2 & 4). Where 
we start to differ appears to be in working out how 
af7?(Ap 4- (5Ap)](Ap -I- dAe) is related to a[fj(Ae)]Ae (see eq. 
[9] of [13 )■ We believe the correct relation (to first order) 
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is 


+ l5Ae)](Ae + 5Xe) (C24) 

= a[77(Ae + l5Ae) + l5?7(Ae)]Ae(l + 6Xe/Xe) 

= a[f7(Ae)]Ae(l+ ^ [5r?(Ae) - <5Ae] + ^) 


The combination 6ri{Xe) — SXe is equivalent to r]{Xe + 
SXe) — fj{Xe), which was referred to as Sr] in Ea , let us 
denote this as SrjPB- Eq. (9) of da can be written as 


®[^(Ae + <^Ae)](Ae + I^Ag) 

V 


(C25) 


= a 


[??(Ae)]Ae( 


1 + 


SrjPB + 


SzpB \ 


where SzpB was defined in eq. (6) of ^3 SzpB = 
X^{Xe + SXe) - X^(Xe) = X^(Xe + SXe) “ X^ (Xe) + Sx^(Xe). 
Using eq. and ra . we can see that SzpB = 

SXe + Sx^(Xe) and ^^(Ae) = -2 dX'(Xe-X')S>,r,- The 

issue is that da appears to have wrongly used SzpB in¬ 
stead of SXe as it should be (compare eq. 1(1241 and 
1(1251 1. and the two differ by Sx^{Xe), which leads to ex¬ 
actly the difference of 2 dX'(Xe — X')/Xe4>,ri in <5di- 
The source of the error in da seems to be confusion in 
the definition of what they refer to as Ve- According to 
an earlier paper by the same authors [T^ . Ve should be 
equated with r(Ae) where r is the comoving radial dis¬ 
tance in an unperturbed, rather than a perturbed, uni¬ 
verse (see eq. [6] and [29] of 0). (Note that our A 
convention is opposite in sign from that in |l(i|: we have 
X® = An* while they effectively have x* = —An*, setting 
their Aq = 0.) 

Eq. EHl) has quite a number of terms. They can 
be loosely divided into four categories: peculiar motion 
(first line), gravitational leasing (second line), gravita¬ 
tional redshift (third line) and integrated Sachs-Wolfe 
(fourth and fifth lines). This division is loose in a number 
of ways. Eor instance, some of what we call the integrated 
Sachs-Wolfe terms (those that involve integrals over time 
derivatives of (j)) can actually be physically ascribed to 
gravitational leasing. We call the second line the gravi¬ 
tational leasing term mainly because the literature often 
uses it as the only contribution to the leasing conver¬ 
gence. The first two terms of the third line, what we call 
gravitational redshift terms, obviously mirror those for 
peculiar motion i.e. Doppler shift, but the last term of 
the third line has no Doppler analog. 

Eortunately, for most applications we can confine our 
attention to the peculiar motion and leasing terms (first 
two lines). This is because we are generally interested 
in fluctuations on scales smaller than the horizon - the 
high redshift SN surveys generally cover a small fraction 
of the sky while the low redshift surveys, even though 
they cover a significant fraction of the sky, do not ex¬ 
tend out to a sufficient depth to be sensitive to horizon 
scale fluctuations. In other words, H/k is a small number 
where H is the Hubble constant, and k is the wavenum¬ 
ber of interest (strictly speaking, the correct quantity to 


look at is aH/k, but here a ~ 1). For instance, the term 
involving (j)" is obviously smaller than the lensing term 
involving V^(/> i.e. (j)" ^ (H/k)^ 1. The rate of 

change of (/>, if non-zero, should be of the order of the 
inverse Hubble time ^ H. Similarly, terms like J dx4>'i 
f dx4>'{Xe — x)IXe, or f (pdx/Xe sre at most of the order 
of (^, which in turn is smaller than v. The latter holds 
because v ~ (k/H)(j) (for a ~ 1), from gravitational in¬ 
stability. 

This is why for the purpose of this paper, we can con¬ 
sider only the peculiar motion and lensing contributions 
to the luminosity distance fluctuation. This justifies the 
use of eq. 


APPENDIX D: THE MAGNITUDE 
COVARIANCE MATRIX II - EXPLICIT 
EXPRESSIONS IN TERMS OF THE MASS 
POWER SPECTRUM 

Our main task here is to derive explicit expressions 
for the lensing and velocity contributions to the (binned) 
magnitude covariance matrix, including both the Poisso- 
nian and non-Poissonian terms i.e. eq. (Cni), 1201), I2U, 
(I22II and At the end of this Appendix, we will also 
give expressions for the velocity window function that are 
useful for more general survey geometries. 

Let us work in the plane-parallel (small angle) approx¬ 
imation first. We will generalize some of the relevant 
terms to large angles later. 

The starting point is the luminosity distance fluctua¬ 
tion SdL given in eq. (111^1 . We substitute this into eq. 
(0 and O to compute ((Tf°*®®')^ and C^. Let us first 
study the lensing contributions. The expressions are, up 
to normalization factors, identical to those worked out 
in the lensing community for convergence variance (e.g. 

dS ), and so our derivation here is very brief. From 
eq. m, using Limber’s approximation and Poisson’s 
equation, it can be shown that 


{SddZi,9)SdAZj,9%ens = (3H2fi^/2)2 (Dl) 
rmin.te,xH dy (x, - x)y (y, - x)x 

Jo o? Xi Xj 

where Xi and Xj are the radial comoving distances (eq. 
I(122p to redshifts Zi and zj respectively, and P{k,a) is 
the mass power spectrum at a scale factor a correspond¬ 
ing to when the photon is at distance y. Here is the 
projection of k onto the plane perpendicular to the line 
of sight. 

It is simple to set 9 = 9' and i = j in the above 
expression and show from eq. that 

2 

{iH^n^/2)^ (D2) 


/ Poiss., lens \2 
Wi ) = 


InlO 












23 


^ r {xi - x)x 

Jo a? [ 


(fk± 


P{k = k±,a) 


provided the redshift bin width is not too large. This 
reproduces the lensing part of the Poissonian variance in 

eq. ((201) • 

We can similarly substitute eq. EU into eq. 0 to 
obtain the lensing term in the non-Poissonian correlation 
matrix in eq. In doing so, it is useful to remember 

that eq. m is a slow function of Zi and Zj, and the 
averages over redshifts within the respective bins in eq. 
o can safely be replaced by the middle values, provided 
the redshift bins are not too large. It is also useful to 
note that for a survey that spans a (not too large) circle 
on the sky with angular radius (area Ai): 


(Pe 

X' 


-ikxxS 


(D3) 
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Jo 
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dxX Inlx) — _^ ^ ^_ L 

axXJo[X) - 


The velocity analogs of the above can be derived by 
using mass conservation (to linear order): 

5' = -V • V (D4) 

where 6 is the density fluctuation Sp/p. This tells 
us that in linear theory and Fourier space: (k) = 

i{D'/D){k^/k^)S(k), where D is the linear growth factor 
and D' is its derivative with respect to conformal time. 
Therefore, we have 


(u3(xi)n3(x2)) = 

f d?k {k^Y 


(D5) 


(27r)3 k'^ 




where we have adopted the z (or third) direction as the 
line of sight direction. The derivatives of the growth fac¬ 
tor D'^ and D 2 are evaluated at redshifts corresponding 
to the positions xi and X 2 . 

The velocity two-point correlation above, together 
with the relation between v and in eq. m, and 
the relation between and SdL in eq. ®, gives 

us the velocity contribution to the Poissonian magnitude 
variance, if we set Xi = X 2 : 


/ Poiss., vel .\2 

) = 

f dJk {kzY 

J (27r)3 k* 


In 10 


(1--X 


? (D6) 


{D[YP{k,a=l) 


which is consistent with eq. ( 1201 ). Note that as in the case 
of lensing, we assume that the redshift bin is sufficiently 
narrow that factors like a, x, D and so on do not vary 
much across the bin, and so one can use their values at 


the center of the bin (i.e. at z = Zi) to substitute for what 
should strictly speaking be bin-averages. In deriving the 
above, we have ignored the Uq term in the luminosity dis¬ 
tance fluctuation (eq. mi): 5di^ = Ug —(a/a'/x)e(Vg—Ug) 
(ignoring the lensing term for now). This is justified to 
the extent that our peculiar motion is fairly well con¬ 
strained by the microwave background dipole, and so one 
could eliminate its influence on luminosity distance mea¬ 
surements. The other extreme would be to allow Ug to 
be a stochastic variable just like in which case the ve¬ 
locity contributions to the magnitude variance would be 
even higher than our estimate. There is one more sub¬ 
tlety: even in the case where Vq is known precisely, our 
computation of {a^ o'ss.. ve .^2 strictly speaking only ap¬ 
proximate. We have essentially computed [5/In 10]^ [1 — 
{a/a'/x)eY{{Ue Y)j whereas the correct thing to do is to 
compute [5/ In 10]2[I - {a/a'/x)ef[{Yl?H) “ ko)^]> 
because there is some correlation between Ug and v/. 
((uglug) is the expectation value of Ug given Ug). If 
the emitter and the observer are sufficiently far apart, 
((Ug)^|uo) ~ {{VeY)- The term (r'glug)^ for Gaussian 
random fluctuations is equal to [(r’er'g)/((r'g)^)]^(r’g)^. 
Approximating (ug)^ by ((ug)^), one can see that our 
computation of j 2 igg fairly accurate if 

(UgUQ)^/[((uQ)^)((z)g)^)] <C 1. We have checked that in 
all cases we have considered in this paper, the corrections 
due to the correlation between Vg and Vg are negligible. 
Similar, but slightly modified, arguments apply to the 
non-Poissonian velocity terms (e.g. the monopole and 
so on) considered below. The relevant quantity in that 
case is ([v(xi) • Xi][vo • Xi])([v(xj) • Xj][vo • Xj])/[([v(xi) • 
Xi][v(xj)-Xj])(vo-vo)], where Xi and xj are to be averaged 
over redshift bins i and j respectively; we have checked 
explicity that this is also ^ 1 for all cases of interest in 
this paper. 

The velocity contribution to the non-Poissonian corre¬ 
lation matrix G™'' can be worked out in an analogous 
manner, and one arrives at eq. 122 ). As before, we have 
made use of the fact that quantities like a, D, x 
so on vary slowly across the redshift bin (but not quan¬ 
tities such as The window function given in eq. 

(22 is for volumes that are top-hat in the z-direction and 
circular in the transverse direction. Note that if z Y 
strictly speaking the window function should have 

an imaginary piece proportional to z sin [kz {Xi~Xj)]i but 
such a piece doesn’t contribute to the Fourier integral 
because it is odd in k^- 

The reader might wonder whether there should also be 
a cross-term that involves the product of velocity-induced 
and lensing-induced fluctuations. Such a term is in prin¬ 
ciple possible, but the lensing projection sends kz to zero 
(Limber approximation) while the velocity contribution 
is proportional to fcj, making the cross-term identically 
zero. (A similar statement holds even in the large angle 
case studied below: essentially, the cross velocity-lensing 
term would involve an integral of je{kx) over a large dis¬ 
tance Xj making it very small.) 

Finally, for realistic applications, we need to generalize 
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the above discussion for velocities to large angles. This is 
because low redshift SN surveys, where peculiar motion 
(but not lensing) is important, typically cover a signifi¬ 
cant fraction of the sky. A more general version of eq. 

m is 


([v(xi) •Xi][v(x2) •X2]) = 

f (Pk (k • Xi)(k • X 2 ) 


(27r)3 


fc4 


(D7) 

D[D'^P{k,a = 


where xi and X 2 are unit vectors pointing towards the 
positions xi and X 2 . It is straightforward to show that 
eq. follows from the above expression plus eq. HU) 
and eq. 0 . It takes a little more work to derive the 
velocity window function given in eq. 12.‘lll for a 

circularly symmetric survey. The definition given in eq. 
(I22II can be written as: 
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1 
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x]AxjAilj 
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X^dxdne-^'^-^k ■ x 


x'^dx'dnp^-^k ■ x' 


where Ay^, Axj are the radial widths and Afli and Aflj 
are the solid angles. 

A useful expansion of the plane wave is: 


g = ^47r(-f)^j^(fcx)r<>„(k)F/„(x) (D9) 


patch retains azimuthal symmetry (i.e. each patch spans 
from (j) = 0 to (j) = 27r; we assume there is no danger of 
confusing the azimuthal angular coordinate (j) with the 
metric fluctuation (j>). This is precisely the kind of geom¬ 
etry we have adopted for the SNfactory, which consists of 
two patches, one centered at the north pole and the other 
at the south pole. For a survey like this, generalizing eq. 
a little bit, the velocity window function should be: 


w:;^-ik) = Y,m + i) 


(D12) 
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where the symbol denotes integration over the patches 
that belong to the Ah redshift bin. For instance, if one 
has a survey that covers all sky aside from a galactic cut 
of ±30°, /. should represent an integration of 9 from 0 to 
7r/3, and from 27r/3 to tt. Eq. 0121), in place of eq. 12211, 
is what we use to obtain predictions for the SNfactory. 

Generalizing further, suppose one has a survey where 
for each angle 9, the azimuthal angle (p spans <pmin.{d) 
to 4>-ma.iL.(d)■ We will further give pram, and ^imax. super¬ 
scripts i or j to denote the fact that these ranges can 
even vary depending on the redshift bin under considera¬ 
tion. It can be shown that the velocity window function 
is given by: 


where ji, is the spherical Bessel function and the are 
spherical harmonics. From this, one can see that 

• X = ^47r(-i)^-^j;(A:x)yf„(k)Y/„(x)(D10) 

where j[ is the derivative of with respect to its argu¬ 
ment. 

Looking at eq. 1221i, one can see that eventually we in¬ 
tegrate kF™* (k) over k with other functions of k (but not 

of k). In other words, we can replace IF™'- (k) by its aver¬ 
age over the solid angle of the wave vector. One can then 
take advantage of the fact that f dOfcy^’(^(k)IY/m'(k) = 
5wdrum'- The integration over the solid angle dVL of the 
survey can be done quite easily: 

J dflYe^ix) = J dnYgo{x.)Smo (DU) 

/ ^max. 

r' 

= 2Try^^ - J d9sm9Pe{cos9)SmO 

The solid angle AXli equals 27r(l — cosfl™'^’^ ). Putting 
all these together yields eq. 1221 ) for the velocity window 
function kF™°, which is valid for a contiguous, circularly 
symmetric survey. 

A slightly more general situation is one where the sur¬ 
vey is divided into several different patches, but each 


VF™' (fc) = ^(2£± l)Aflr^AfIj 


(D13) 
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An^ = [d9 sine 
Jj 


and P™ denotes the associated Legendre polynomial. 
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